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1.0  INTRODUCTION 

This  report  describes  progress  during  the  third  year  of  a  research 
program  on  the  Fundamental  Processes  in  Partially  Ionized  Plasmas  con¬ 
ducted  in  the  High  Temperature  Gasdynamics  Laboratory  at  Stanford  Uni¬ 
versity.  This  research  is  supported  by  a  grant  from  the  Air  Force  Office 
of  Scientific  Research  (AFOSR-83-0108)  and  is  conducted  under  the  direc¬ 
tion  of  Professors  Charles  H.  Kruger,  Morton  Mitchner,  and  Sidney  Self. 
Three  Ph.D.  candidates  are  currently  conducting  their  doctoral  research 
under  this  program. 

Several  space  power  and  propulsion  systems  of  potential  long-range 
interest  to  the  Air  Force  involve  partially  ionized  plasmas.  Such 
systems  include  MPD  thrusters,  both  open  and  closed  cycle  MHD  power,  and 
thermionic  energy  conversion.  Although  the  specific  configurations,  the 
exact  operating  conditions,  and  which  of  the  competing  systems  will 
prove  to  be  most  useful  in  the  long  term  remain  to  be  established,  it  is 
important  at  this  time  to  provide  a  broad  fundamental  research  base  in 
support  of  development  activity.  In  particular,  there  are  a  number  of 
key  issues  regarding  the  properties  and  discharge  behavior  of  partially 
ionized  plasmas  and  the  interaction  of  discharges  with  fluid  dynamics 
that  need  to  be  understood  before  the  potential  and  limitations  of 
competing  systems  can  be  fully  evaluated.  In  addition,  it  is  important 
that  outstanding  young  applied  scientists  be  educated  in  these  areas. 

The  present  research  on  partially  ionized  plasmas  is  discipline 
rather  than  device  oriented  and  is  currently  focused  on  three  major 
areas : 

1.  Recombination  in  molecular  plasmas 

2.  Discharge  effects:  plasma  electrode  interaction 

3.  Interaction  of  discharges  and  fluid  dynamics 


In  addition,  each  of  these  areas  involves  the  development  and  applica¬ 
tion  of  modern  plasma  diagnostic  techniques. 


These  areas  are  overlapping  and  mutually  supportive.  Thus,  under¬ 
standing  of  plasma  properties  is  important  to  the  study  of  discharges 
and  their  interaction  with  fluid  dynamics.  In  the  same  spirit,  we  are 
interested  in  the  development  of  plasma  diagnostics  so  that  they  can  be 
applied  in  our  research  objectives.  In  each  area,  the  research  is  pri¬ 
marily  experimental  in  nature,  with  supporting  theoretical  studies  for 
the  planning  of  the  research  and  interpretation  of  the  data. 

Progress  during  the  third  year  in  each  of  the  three  research  areas 
is  described  by  means  of  Project  Summaries  in  Section  2.0. 

Publications  and  Presentations  resulting  from  this  work  are  cited  in 
Section  4.0,  and  Section  5.0  lists  the  personnel  who  have  contributed  to 
this  report. 


2. 1  Recombination  in  Molecular  Plasmas 
Introduction 

In  applications  of  interest  to  the  Air  Force  (e.g.  MHD  generators 
and  space  propulsion)  where  a  high  pressure  (i.e.  of  the  order  of  an 
atmosphere)  partially-ionized  plasma  flows  past  a  cooled  surface,  recom¬ 
bination  of  electrons  occurs  with  positive  alkali-metal  ions  (such  as 
Cs+)  under  conditions  for  which  there  have  been  few,  if  any,  reported 
measurements  of  the  recombination  rate.  Typically,  these  conditions  are 
for  relatively  low  electron  number  densities  ng  where  the  three-body 
recombination  process 

e  +  Cs+  +  e  -►  Cs  +  e  ,  (l.l) 

1  2 

usually  dominant  at  higher  electron  number  densities  (e.g.  ng  £  10 
cm"\  at  T  2  *1000  K  ) ,  no  longer  plays  a  major  role. 

To  obtain  an  improved  understanding  of  recombination  processes  under 
such  conditions,  we  have  proposed  to  design  and  perform  experiments  in 
which  a  highly  nonequilibrium  plasma  is  produced  in  a  mixture  containing 
cesium  vapor  and  nitrogen  gas.  From  an  engineering  point-of-view,  it  is 
the  actual  value  of  the  recombination  rate  which  is  of  primary  impor¬ 
tance,  independent  of  what  the  underlying  kinetic  mechanisms  may  be. 
However  a  knowledge  of  the  mechanisms  is  very  useful  in  that  it  enables 
one  to  extrapolate  and  predict  values  of  the  recombination  rate  for 
conditions  that  go  beyond  those  for  which  measurements  exist. 

In  addition  to  measurements  of  the  recombination  rates,  another 
objective  of  our  proposed  work  is  therefore  to  identify  the  operative 
recombination  mechanisms  and  to  make  comparisons  of  our  measurements 
with  existing  theories.  In  previous  descriptions  of  our  proposed  work, 
we  have  discussed  the  three-body  recombination  process 

e  +  Cs+  +  N2  -v  Cs  +  N2  (1.2) 

as  one  that  would  occur  under  conditions  of  interest  (e.g.  a  MHD  bound¬ 
ary  layer).  The  recombination  coefficient  for  this  reaction  has  been 


calculated  by  Bates  et  al.  [1.3],  but  no  experiments  have  been  per¬ 
formed  to  date  that  would  enable  a  comparison  to  be  made  with  theory. 


One  experimental  arrangement  for  producing  the  nonequilibrium  condi¬ 
tions  which  are  necessary  to  study  recombination  is  to  employ  a  super¬ 
sonic  nozzle  to  rapidly  expand  a  flowing  high  temperature  plasma,  con¬ 
sisting  of  N2  and  thermally  ionized  Cs.  Measurements  of  the  electron 
number  density  before  and  after  the  expansion  would  be  compared  with 
calculated  values  obtained  from  the  electron  continuity  equation.  Since 
the  calculated  values  of  ne  depend  on  the  recombination  coefficient 

aM  (T),  such  an  experiment  would  provide  a  means  for  inferring  the 
"2 

recombination  rate.  (Here  a  is  defined  for  recombining  conditions  by  the 
relation  dng/dt  =  -  ang  n^,  where  n^  denotes  the  number  density  of  the 
ions  with  which  the  electrons  are  recombining.) 

To  analyze  more  accurately  the  feasibility  of  such  an  experiment,  a 
mathematical  model  of  the  flow  was  formulated,  and  solved  numerically. 
The  results  showed  that  in  order  for  the  value  of  ng  at  the  exit  of  the 
nozzle  to  depart  significantly  from  its  frozen  value,  the  upstream  value 
of  ng  would  need  to  be  sufficiently  large  that  the  process  (1.1)  would 
contribute  appreciably  to  the  over-all  recombination  rate.  Thus  the 
accuracy  with  which  could  be  determined  would  be  reduced.  This  kind 
of  an  experiment  also  has  the  disadvantage  that  the  recombination  would 
be  occurring  under  non-isothermal  conditions,  thereby  further  compli¬ 
cating  the  reduction  of  the  data. 

Because  of  these  limitations  (as  well  as  others)  effort  has  been 
directed  toward  analyzing  alternative  experimental  arrangements  for 
producing  the  required  nonequilibrium  conditions.  Among  the  alternatives 
considered,  the  approach  that  appeared  most  attractive  was  one  in  which 
the  ionizational  nonequilibriura  would  be  achieved  by  pulsed  photoioniza¬ 
tion  of  Cs  vapor  confined  with  N2  in  a  heated  test  cell.  The  recombina¬ 
tion  coefficient  aN^(T)  would  be  determined  by  measuring  the  decay  of  ng 
as  a  function  of  time,  following  the  completion  of  the  photoionization 
pulse.  In  a  photoionization  experiment  the  recombination  is  non-steady 
in  time  and  uniform  in  space,  in  contrast  to  a  rapid  expansion  experi¬ 
ment  which  is  steady  in  time  and  nonuniforra  in  space. 


In  the  remainder  of  this  section,  we  outline  first  our  research 
objectives  and  then  discuss  the  present  status  of  the  research.  During 
the  past  year  a  theoretical  analysis  directed  toward  showing  in  more 
detail  the  feasibility  of  a  pulsed  photoionization  experiment  was  com¬ 
pleted,  and  substantial  progress  was  made  toward  completion  of  the 
experimental  facility.  The  main  result  of  our  present  feasibility  study 
is  to  show  that  for  mixtures  of  interest,  below  a  temperature  of  about 
800K,  recombination  involves  several  ion-molecule  reactions  proceeding 
simultaneously.  (See  the  Feb.  1,  1985  to  Jan.  31,  1985  Annual  Scientific 
Report  for  a  description  of  past  feasibility  calculations.)  However  for 
temperatures  in  excess  of  about  1000K,  electron  recombination  is  domi¬ 
nated  by  the  three-body  process  (1.2).  Our  recombination  rate  measure¬ 
ments  at  higher  temperatures  should  therefore  enable  us  to  make  compar¬ 
isons  with  the  Bates  theory.  Also  during  the  past  year  much  of  the 
experimental  facility  has  been  fabricated.  Items  completed  include  the 
reaction  cell,  the  heating  and  vacuum  system,  and  a  prototype  flash  lamp 
power  supply. 

Research  Objectives 

1.  To  design  and  construct  an  experimental  system  that  can  produce  a 
recombining  alkali  metal  plasma  under  conditions  where  the  molecular 
three-body  process  dominates. 

2.  To  examine,  in  particular,  the  recombination  of  cesium  ions  in  the 
presence  of  nitrogen. 

3.  To  measure  the  desired  recombination  rate  as  a  function  of  the 
background  gas  temperature. 

4.  To  measure  the  recombination  rate  as  a  function  of  various  values  of 
the  partial  pressures  of  Cs  and  N2 • 
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5.  To  examine  the  recombination  rate  of  interest  in  a  mixture  of  molec¬ 


ular  species. 

b.  To  design  and  construct  a  suitable  pulse  forming  circuit  for  a  high 
intensity  flashlamp  and  to  explore  the  use  of  pulsed  photoionization 
as  a  means  for  producing  the  desired  plasma  conditions. 

7.  To  determine  the  electron  number  density  as  a  function  of  time  in  a 
plasma  afterglow,  and  to  deduce  the  recombination  rate  from  this 
measurement . 

8.  To  develop  a  four  pin  probe  for  measuring  the  time  resolved  plasma 
electrical  conductivity  as  a  diagnostic  for  determining  the  electron 
concentration. 

9.  To  develop  a  method  employing  an  equilibrium  plasma  for  calibrating 
the  electron  conductivity  probe. 

10.  To  examine  the  possible  use  of  optical  diagnostic  methods  such  as 
measurement  of  the  emission  of  free-bound  continuum  radiation,  as  a 
means  for  obtaining  a  check  on  the  data  via  an  independent,  redun¬ 
dant  proceedure. 

11.  To  assess  the  possible  roles  of  the  presence  of  the  dimer  ions  Cs9+ 
and  ion  clusters  such  as  Cs4"*!^  and  Cs^f^O  (the  water  vapor  will  be 
present  in  small  quantities  as  an  impurity)  in  contributing  to  the 
measured  recombination  rates. 

Status  of  Research 

During  the  period  February  1985  through  January  198b  research 
progress  consisted  primarily  of  the  development  of  a  detailed  computer 
model  of  the  reaction  kinetics  for  the  pulsed  photoionization 
experiment,  and  of  the  construction  and  assembly  of  major  elements  of 
the  experimental  facility.  The  model  includes  important  ion-molecule 


chemistry,  and  it  demonstrates  the  validity  of  electron  loss  measure¬ 
ments  as  a  means  for  determining  the  rate  of  reaction  (1.2).  4  reaction 
cell,  three  ^onal  ovens,  a  temperature  measurement  system,  a  temperature 
control  system,  and  a  vacuum  system  were  all  completed  during  the  past 
year.  Additionally,  a  flashlamp  power  supply  was  prototyped. 

Reaction  Kinetics  Model 


In  previous  discussions  we  have  demonstrated  that  experimental 
conditions  can  be  chosen  such  that  recombination  by  process  (1,2)  would 
dominate  that  by  process  (1.1)  In  addition  to  process  (1.2),  a  loss  of 
electrons  under  conditions  of  interest  could  also  occur  through  a 
sequence  of  reactions  involving  fast  formation  of  the  ion  cluster 
Cs+*N2,  fast  2-body  switching  to  form  Cs^,  followed  by  dissociative 
recombination  of  the  cesium  dimer  ion.  These  processes  can  be  written  as 
follows : 


Cs+  +  N2  +  N2  *  Cs+*  N2  +  N2  (1.3a) 

Cs+*N2  +  Cs  *•  Cs*  +  N2  (1.4a) 

e  +  Cs*  >  Cs  +  Cs  (1.5) 

Reactions  analogous  to  (1.3a)  and  (1.4a)  can  occur  when  water  vapor 
impurities  are  present  and  H2O  replaces  ^2  as  a  clustering  partner. 
These  reactions  are  written: 

Cs*  +  H20  +  N2  +  Cs*  »20  +  M2  (1.3b) 

Cs+*H20  +  Cs  *•  Cs*  +  H70  (1.4b) 

Other  possible  mechanisms  of  electron  loss  are  electron  dissociative 
recombination  with  Cs+*  N2,  electron  attachment,  and  ion-ion  recombina¬ 
tion  as  represented  respectively  by  the  reactions 


e  +  Cs+*N7  *  Cs  +  N7 


(1.6) 


e  +  Cs  +  N,  ♦  Cs  +  N„ 


(1.7) 


Cs  +  Cs  +  2Cs 


(1.8) 


In  what  follows,  we  present  a  detailed  examination  of  the  overall 
recombination  process  that  includes  both  the  ion-molecule  reaction 
channel  (1.3),  (1.4),  and  (1.5)  as  well  as  the  three-body  processes  (1.1 
and  1.2).  Using  estimates  of  the  recombination  rate  coefficients  and 
equilibrium  constants  (to  calculate  the  backward  rates)  constructed  from 
available  information  in  the  literature  (see  Table  1.1  for  references), 
we  show  that  experimental  conditions  at  higher  temperatures  can  be 
chosen  to  maintain  process  (1.2)  as  the  dominant  electron  loss  mech¬ 
anism,  despite  the  fast  reactions  (1.3)  and  (1.4).  Later  we  show  that 
inclusion  of  reactions  (1.6) — ( 1.8)  does  not  affect  the  dominance  of 
reaction  (1.2)  in  electron  loss  frequency.  In  addition  we  also  show  that 
by  using  a  liquid  nitrogen  cold  trap  and  because  of  the  gettering  action 
of  Cs,  as  well  as  from  considerations  of  the  reaction  kinetics,  the 
deleterious  effects  of  a  possible  water  impurity  are  not  significant. 

Partial  Equilibrum 

Although  reactions  (1.3)  and  (1.4)  proceed  in  a  forward  direction 
rapidly  relative  to  the  three  other  reactions  considered,  the  backward 
rates  for  these  reactions  are  also  large.  Reactions  (1.3)  and  (1.4)  are 
able  to  quickly  reach  a  partial  equilibrium  condition  where  the  forward 
and  backward  rate  are  equal  in  magnitude.  This  partial  equilibrium  is 
maintained  while  the  recombination  reactions  slowly  change  the  ion  and 
electron  concentrations. 

Table  1.1  lists  the  preceeding  ten  chemical  reactions.  Also  listed 
are  the  values  of  the  corresponding  forward  and  backward  rate  constants 
and  the  ratios  of  these  two  constants  (i.e.,  the  equilibrium  constants), 
used  in  all  subsequent  calculations.  The  constants  were  obtained  using 
the  references  shown.  The  photoionization  of  Cs  vapor  by  ultraviolet 
radiation  from  a  flashlarap  can  be  represented  by  the  reaction 


V  + 

Cs  +  hv  -*■  Cs  +  e  . 
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Reaction 


+ 

3  +  e 

+• 

k 


k 

+  > 


Table  1.1  Reactioa  Rate  and  Equilibrium  Coastauts 


Constants 


References 


k^i.ogxio"8/?4*5 

Cs  +  e 

k_1-2.64X107/T3*e 

■1 

kl/k_l=4.13X10"16 

2 

k2=,6 .58X10~29*exp 

’  Cs  +  N2 

k_2=l . 59X10” 13*T3 

-2 

k2/k-2=4- 13x10“16 

:3a  + 

Cs  •  N2  +  N2 

k3a=6  .42X10“^  cm 
k-3a=7*81xl°18*T~ 

-3a 

k3a'/k-3a=8*22xl°" 

:3b  + 

Cs  •  H20+N2 

k3b=»6.42X10~10  cm 
k_3b=2.44X1018/T*i 

-3b 

k3b/k-3b=2,63xl0_ 

1.9, 1.1 


1.15,1.16,1.18 

1.10,1.11,1.12 


1.16,1.8,1.18 


+  >4a  + 

Cs  •  +  Cs  *  Cs^  +  t*2 

k-4a 


k4a=10 


=in-9  cm3/ j 


+  +4b  + 

Cs  •  H20  +  Cs  *  Cs2+  H20 

k-4b 


k-4a=exP(_1 *08 1X10  /T~23* 97+8 •303X10~4*t"  1*17,1.19,1.20 
-1.496X10-7*T2)  cm3/ s 
k4a/k-4a  =  exp(  1  .08 1 X104/T  +  3.244- 
-8.303X10-4*T+1 .496X10-7*T2) 


k4b=*10  9  cm3/s 

k_4b=exp(-7.134X103/T-23.39+8.525X10~A*T-  1.17,1.19,1.2'' 
-1.559X10_7*T2)  cm3/ s 
k4b/k-4b  =  exp(7.134X103/T+2.668- 
-8.525X10-4*T+1.559X10-7*T2) 


Table  l.i  Reaction  Rate  and  Equilibrium  Constants  (continued) 


Reaction 


Constants 


References 


k 

s 

k5-10  *(300/T)  2  cm  /s 

e  + 

+ 

Cs~  Cs  +  Cs 

k_5=exp(-3.094X104/T+2.5*ln(T)-41 .57-  1 

.21,1.22 

k-s 

-8.303X10"A*t+1 .496X10“7*T2)  cm3/s 

k5/k_5-exp(3.094X104/T--3*ln(T) 

+26.00+8.303X10~4*T  -  1 .496X10~7*T2) 

k 

k8» 5.0X1 0_7*( 300. 0/T)  cm3/ s 

e  + 

+ 

Cs  *N2  Cs  +  No 

k_6-5.2X10“16*T1 ,5*exp(-41754./T)  cm3/s  1 

.22 

k-6 

k6/k-6”5,02X10l2*T~3*exP(41754*/T) 

i?  . 

Cs  +  N2  ^  Cs  +  Nj 

k7=8.2X10  *4  cm3/s 

e  + 

k_7-1.98X102*T1'5*exp(-5452.1/T)  1/s  1 

.31 

k-7 

k/7/k_7-4.14X10-i6*T~1-5*exp(5452.1/T)  cm3 

k8 

k8-2.0X10_7*300/T  cm3/ s 

Cs~ 

+  Cs  ^  Cs  +  Cs 

k_8-2.0X10"7*300/T*exp(-39673/T)  cm3/s  1 

.11,1.32,1.33 

k-8 

ka/k  Q=*exp(  39673/T) 

^  o 

The  photon  flux  is  modeled  as  a  rectangular  function  in  time,  and  the  Cs 

vapor  can  be  shown  to  be  optically  thin.  Under  these  conditions  the  rate 

of  increase  of  the  electron  number  density  can  be  represented  by 

d[e]/dt  ■  [Cs]ky  where  kv  is  the  integral  over  wavelength  of  the  product 

of  the  photon  flux  and  the  total  photoionization  cross  section.  The 

value  for  kv  used  in  these  calculations  is  1.0  sec  This  value  corre- 

20 

spends  to  a  photon  flux  (wavelengths  <  318  nm)  of  about  10 

2 

photons/cm  *sec,  which  is  an  estimate  of  the  flux  attainable  with  our 
pulsed  short  arc  flashlamp  system. 

Shown  in  Figs.  1.1a,  1.2a,  and  1.3a  are  the  time  histories  of  the 
number  densities  of  the  species  of  interest  for  three  different  test 
cell  temperatures.  The  Cs+*H20  species  occurs  in  concentrations  which 
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Fig.  1.2a  Species  concentrations  from  numerical  solution  of  kinetics 
equations  for  reactions  ( 1 )  — ( 5)  .  T>1200  X;  [Cs]=4.5X10*4  cm  3;  [N2]»10*^ 
cm--*;  [H2o]«1013  cm-3.  Two  distinct  regimes  are  shown:  t<0,  complete 
equilibrium;  t>l  ms,  plasma  relaxation.  Photoionization  occurs  between 
these  regimes. 
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Fig.  1.2b  Electron  loss  frequencies  from  miner leal  solution  of  kinetics 
equations  for  reactions  (l)-(5).  T»1200  K;  [ Cs ] *4 .5X10  * ^  cm  3; 
[  V2  ]  —  1 0  ^  ^ ;  (H20]**10*3  cm-3.  l/r^(e)  i  (net  electron  loss  rate)/[e].  For 
ambipolar  diffusion,  l/x^lff(e)  £  10  sec-*. 
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Fig.  1.3a  Species  concentrations  from  numerical  solution  of  kinetics 
equations  for  reactions  ( 1)  —  (5)  •  T“1500  K;  [ Cs ) “4 . 5X101 ^  cm  3;  (N2)“1019 
cm--*;  [H2O]=1013  cm-3.  Two  distinct  regimes  are  shown:  t<0,  complete 

equilibrium;  t>l  ms,  plasma  relaxation.  Photoionization  occurs  between 


these  regimes. 
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Fig.  1.3b  Electron  loss  frequencies  from  numerical  solution  of  kinetics 
equations  for  reactions  (l)-(5).  T^ISOO  K;  [  Cs  ]  “4.5X101 14  cm  3; 
[N2]-1019;  [ H20 ] = 1 0 1 3  cm-3.  l/t^(e)  =  (net  electron  loss  rate)/[e).  For 
ambipolar  diffusion,  1/Ttlff(e)  £  10  sec-1. 


are  off  the  lower  end  of  the  plotted  scale.  This  species  is  treated 
explicitly  in  the  section  on  impurities.  (Note  that  reactions  (1.6)  - 
(1.8)  are  not  modeled  in  the  results  presented  in  these  figures.)  It  is 
assumed  that  the  flash  lamp  is  energized  at  t=0  and  it  is  extinguished 
at  t=l  ms;  the  dynamic  relaxation  of  the  nonequilibrium  plasma  occurs 
for  times  t>l  ms.  These  plots  were  made  by  numerically  solving  the  set 
of  differential  equations  which  describe  the  kinetics  of  Che  species 
involved  in  reactions  (1.1)  through  (1.5).  For  t<0,  the  plasma  is  in  a 
state  of  local  thermodynamic  equilibrium.  The  equilibrium  species 
concentrations  shown  in  the  figures  are  determined  by  the  Cs  vapor 
pressure  (controlled  by  the  reservoir  temperature) ,  by  the  initial  N2 
concentration,  and  by  the  temperature  of  the  reaction  chamber.  Typical 
values  for  these  quantities  were  taken  to  be  [Cs]=*4.5X101,!f  cm-3, 
[N2]*1019  cm-3,  [H20]»1013  cm-3,  and  600  K  <  T  <  1500  K. 

For  a  general  reaction  i  where 

ki 

A  +  B  ♦  C  +  D 

we  can  obtain  order  of  magnitude  estimates  of  the  characteristic  times 
of  significant  decrease  or  increase  of  concentration  of  species  respec¬ 
tively  by,  for  example,  rj[(A)  ~  (kjjB])-1,  or  t£(C)  ~  [C  ]  (ki  [  A]  [3  ]  )-1 . 
In  what  follows,  it  is  shown  how  relations  such  as  these  may  be  employed 
to  illustrate  the  essential  physics  underlying  the  numerical  solutions 
presented  in  Figs.  1.1a,  1.2a,  and  l.3a.  Impurity  reactions  (1.3b)  and 
(1.4b)  have  been  included  in  the  numerical  solutions  shown  in  Figs.  1 .  I  — 
1.3,  hut  an  explanation  of  the  effect  of  including  these  reactions  is 
deferred  until  the  section  on  impurities. 

Using  the  equilibrium  concentrations  for  t<0  shown  in  Fig.  l.2a, 
characteristic  times  for  significant  changes  in  species  concentrations 
are  calculated  in  Table  1.2,  for  a  temperature  of  1200  K.  Considering 
the  Important  reactions  affecting  the  electron  concentration,  the 
ordering  Tpu^ge  =  10-3  sec.  <<  t[<  t^<  demonstrates  that  the  duration 
of  the  photoionization  reaction  (one  millisecond)  is  sufficiently  short 
to  throw  the  plasma  out  of  equilibrium.  The  recombination  reactions 


Table  1.2  Characteristic  Times 
t-0,  T=*1200  K,  Species  concentrations  from  Fig.  1.2a 


Tpulse  duration  =  ^  sec 


T2(e) 

=•  ( [Cs+]  *[N2]  .k2)-1  -  (109«1019-10-28)-1 

=*  1  sec 

T^(e) 

-  ( [Cs+2]  *k5)-1  =■  (3X106*5X10-9)-1  -  102 

sec 

t  J(e) 

■  ( [Cs+] • [e ] *kj )-*  =  ( 109  •  109  •  10-22)_1  =* 

104  sec 

4*<Cs+> 

-  UN2]*k3a)~1  -  (1019*6X10“10)”1  -  lCf9 

sec 

x2(Cs+) 

■  ( [e]  •  [N2]  •k2)~1  =•  ( 109  *10^9  •  10-2®)-1  - 

1  sec 

T4a^Cs+2^  *  [Cs+2]/( [Cs+-N2] *[Cs] *k4a)  » 

=*  3X106/(8X105*5X10U*10-9)  =  10“5  sec 

T^( Cs+2)  ■  ([e]*k5)_1  *  CIO9 • 5X1 0“ 9 )~ 1  ■  1  sec 


(1.1),  (1.2),  and  (1.5)  will  not  affect  the  change  in  ion  and  electron 
concentrations  during  and  immediately  after  the  flash  because  these 
reactions  deplete  charged  particle  concentrations  in  a  characteristic 
time  which  is  much  longer  than  the  pulse  duration. 

The  orderings  T2a(Cs+)<<T2(Cs+)  and  T^a(Cs+2)<<T^(Cs+2)  [note  that 
these  t's  are  calculated  at  equilibrium  where,  for  any  species  s,  tJ^(s) 
■  T^(s)]  indicate  that  reactions  (1.3a)  and  (1.4a)  will  proceed  to 
establish  partial  equilibrium  concentrations  of  Cs+*N2  and  Cs%  while 
the  recombination  reactions  change  the  ion  and  electron  concentrations 
relatively  slowly. 

This  partial  equlibrium  is  embodied  in  the  relations 

k3, 

I—  '  -  — as  - 

[Cs+1[NJ  K 


(1.9) 


Equilibrium  relation  (1.10)  corresponds  to  the  reaction  Cs^  *  Cs  +  Cs+ 
obtained  by  adding  reactions  (1.3a)  and  (1.4a).  Once  reactions  (1.3a) 
and  (l.4a)  are  shown  to  be  fast  enough  to  maintain  partial  equilibrium, 
the  equilibrium  relations  need  not  correspond  directly  to  the  reaction 
channel.  Relations  (1.9)  and  (1.10)  have  been  chosen  to  express  the 
equilibrium  condition  in  order  to  illustrate  that  the  equilibrium  con¬ 
centrations  of  Cs^  and  Cs"^  N2  are  approximately  decoupled. 

Since  N2  and  Cs  are  present  in  much  larger  quantities  than  Os'*" ,  Cs,, 
or  Cs+*N2  (as  can  be  justified  by  the  results  of  the  calculation  being 
described),  the  following  major  species  relations  can  be  written: 

[N_]  •  constant  »  10^  cm  3,  (1.11) 
and  ,  , 

[Cs]  =  constant  ■  4.5  X  101  cm  .  (1.12) 

From  Equations  (1.9)  and  (1.11)  and  Table  1.1  the  ratio  [Cs+» N2 ] / [Cs*]  ■ 
1^2^ ‘k3a/^l-3a  can  be  shown  Co  be  less  than  1*  at  all  temperatures 
greater  than  500  K.  Additionally,  it  has  been  shown  that  the  recombina¬ 
tion  reactions  are  too  slow  to  significantly  affect  charge  concentration 
during  the  flash.  Neglecting  the  Cs+*N2  concentration  and  the  effects  of 
recombination,  charge  conservation  during  the  flash  can  be  written 

[Cs+]  +  [Cs+J  -  [+]  (1.13) 

where  [+]  denotes  the  concentration  of  positive  (or  negative)  charge 
produced  by  the  photoionization  process.  The  charge  produced  by  the 
flash  is 

[  +  ]  =  [Cs]*k6*At  =■  4.5X1014*1*10-3  =  4.5X1011  cm-3  (1.14) 

Equations  (1.10),  (1.12),  (1.13),  and  (1.14)  are  sufficient  to  determine 
equilibrium  concentrations  of  Cs^  and  Cs+;  the  value  for  [Cs+]  can  then 
be  used  with  Equations  (1.9)  and  (1.11)  to  determine  the  Cs  +  *.N2  concen¬ 
tration.  The  species  concentrations  obtained  in  this  way  are  shown  by 


the  curves  plotted  in  Fig.  1.4,  and  compare  well  with  the  results  of  the 
numerical  solution  at  t=l  ms  which  have  been  plotted  as  individual 
points  over  these  curves.  The  results  shown  in  Fig.  1.4  can  be  used  to 
predict  the  relative  rates  of  the  recombination  reactions  which  follow 
the  flash. 

Characteristic  Times  for  Recombination 

Shown  in  Fig.  1.5  are  the  characteristic  reciprocal  times  1/t  for 
significant  electron  recombination  by  each  of  the  processes  (1.1), 
(1.2),  and  (1.5),  at  t=l  ms,  immediately  following  termination  of  the 
flash.  The  quantity  1/t  is  a  measure  of  the  frequency  of  electron  loss 
via  a  given  recombination  reaction.  The  continuous  curves  were  calcu¬ 
lated  using  the  formulas  shown  in  Table  1.2  evaluated  for  the  concentra¬ 
tions  given  in  Fig.  1.4.  At  temperatures  above  1000  K  reaction  (1.2)  is 
clearly  the  dominant  loss  mechanism  immediately  after  the  flash. 

The  numerical  solutions  shown  in  Figs.  1.1b,  1.2b,  and  1.3b  give 
detailed  information  about  electron  loss  mechanisms  for  extended  times 
after  the  flash.  These  plots  correspond  in  temperature  and  time  to  Figs. 
1.1a,  1.2a,  and  l.3a.  The  frequency  of  electron  loss  by  reaction  i, 

_ 1 _  _  (forward  -  backward  reaction  rate) 

(e)  [e] 

is  shown  for  recombination  reactions  (1.1),  (1.2),  and  (1.5)  during 
plasma  relaxation.  The  net  reaction  rate  (forward  -  backward  rate)  must 
be  considered  for  t>l  ms  when  the  reactions  can  approach  equilibrium.  At 
t=l  ms  the  net  rate  and  forward  rate  are  approxiamately  equal.  Generally 
the  electron  concentration  and  the  net  reaction  rate  change  with  time, 
hence  this  frequency  is  usually  time  dependent. 


A  recombination  reaction  is  considered  to  be  the  dominant  electron 
loss  mechanism  when  its  electron  loss  frequency  exceeds  the  frequencies 
of  all  other  competing  loss  mechanisms.  If  measurement  of  electron 
concentration  decay  is  to  be  useful  in  determining  the  rate  of  reaction 
(1.2),  experimental  conditions  must  chosen  so  that  reaction  (1.2)  is  the 
dominant  electron  loss  mechanism.  For  all  temperatures  greater  than 
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FLg.  1.4  Partial  equilibrium  species  concentrations  at  t=l  ms  (immed¬ 
iately  after  photoionization) .  Continuous  lines  correspond  to  analyt¬ 
ical  solution  of  equations  (1.10) — < 1.16);  discrete  points  correspond  to 
numerical  solutions  of  kinetics  equations  as  in  Figs.  1.1-1. 3. 
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Fig.  1.5  Electron  loss  frequencies  at  t=lms  (immediately  after  photo¬ 
ionization).  Continuous  lines  correspond  to  the  formulae  of  Table  1.2 
and  the  concentrations  of  Fig.  1.4.  Discrete  points  correspond  to  numer¬ 
ical  solutions  as  in  Figs.  1.1-1. 3.  For  ambipolar  diffusion, 
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.-1 


1.1-1. 3. 


1000K,  and  all  times  shown,  reaction  (1.2)  dominates  reactions  (1.1)  and 
(1.5)  in  electron  loss.  The  electron  loss  frequency  due  to  ambipolar 
diffusion,  (at  about  1200  K  and  one  atmosphere  of  J^)  is 

approximately  constant  at  10  sec-*.  As  the  electron  concentration  decays 
with  time,  l/T^e),  decreases  because  the  reactions  move  toward  an  equi¬ 
librium  balance  where  the  forward  and  backward  rates  are  equal.  Electron 
concentration  measurements  are  only  useful  while  the  resulting  reduced 
electron  loss  frequency  remains  larger  than  1/t^j^ f (e) ,  hence  for  the 
conditions  of  Figs.  1.1-1. 3  electron  concentration  measurements  should 
be  made  in  approximately  the  first  10  ms  of  recombination.  Electron  loss 
frequencies  at  one  ms  which  are  taken  from  numerical  solutions  similar 
to  those  shown  in  Figs.  1.1b,  1.2b,  and  1.3b  are  plotted  as  discrete 

points  with  the  the  analytical  solution  shown  in  Fig.  1.5.  Both  solu¬ 
tions  agree  that  at  temperatures  above  1000  K,  during  the  first  ten 
milliseconds  after  photoionization,  measurement  of  electron  concentra¬ 
tion  decay  will  provide  direct  information  about  reaction  (1.2). 

Reaction  1.6  competes  directly  with  reaction  1.2  for  electrons,  but 
the  electron  loss  frequency  can  be  expected  to  be  low  due  to  its 
dependence  on  the  low  concentration  of  Cs"1”*^.  Reactions  1.7  and  1.8 
together  constitute  a  channel  for  electron  loss  via  recombination  with 
Cs+.  These  three  reactions  have  been  included  in  an  extended  computer 
model  of  the  kinetics.  The  results  of  the  numerical  calculation  of 
concentrations  and  electron  loss  frequencies  based  on  reactions  (1.1)  - 
(1.8)  are  shown  in  Figs.  1.6a  and  1.6b  (and  can  be  compared  with 
Fig.  1.2).  Figure  1.6b  indicates  that  reactions  (1.6)  and  (1.7)  are 
significant  secondary  electron  loss  mechanisms,  but  the  primary  electron 
loss  mechanism  remains  reaction  (1.2).  Measurement  of  the  rate  of  elec¬ 
tron  concentration  decay  will  be  a  valid  indication  of  the  rate  of 
reaction  (1.2). 

Impurities 

The  species  Cs ,  Cs  +  ,  e,  and  ^  are  necessary  for  the  experiment; 
other  chemical  species  which  may  be  formed  from  these  cannot  be 
avoided.  In  addition,  the  presence  in  the  test  cell  of  impurities  due  to 


Fig.  1.6a  Species  concentrations  from  numerical  solution  of  kinetics 

equations  for  reactions  (l)-(8).  T=1200K.  The  following  initial  species 

1 4  “3  19  —  3 

concentration  remain  roughly  constant  [Csl=4.5xl0  cm  ;  [ N2 ] = 1 0  cm  ; 
[H20=l-0^cm'"^«  Two  distinct  regimes  are  shown:  t<0,  complete  equi¬ 
librium;  t>l  ms,  plasma  relaxation.  Photoionization  occurs  between  these 
regimes . 


Fig.  1.6b  Electron  loss  frequencies  from  numerical  solution  of  kinetics 
equations  for  reactions  ( 1)— (8).  T=1200K.  The  following  initial  species 
concentrations  remain  roughly  constant  [  Cs  J=4 . 5xl0^cm  I.^J^IO^  ; 
[H->0]  =  10I3cm-3.  1/tj  i  (d[e]/dt)i/(e) . 


outgassing  or  impure  reactant  gases  must  be  considered.  Impurity 
concentrations  can  easily  be  kept  small  enough  to  ensure  that  N?  is  the 
dominant  neutral  partner  in  all  reactions,  but  small  amounts  of  radical 
species  could  substantially  affect  ion-molecule  chemistry. 

Water  vapor  is,  of  all  impurities,  most  likely  to  be  both  suffi¬ 
ciently  reactive  and  plentiful  to  have  a  significant  effect.  The 

amount  of  l^O  present  in  the  N2  charge  will  be  minimized  by  passing 
the  gas  through  a  liquid  nitrogen  trap.  The  partial  pressure  of  the 
H2O  will  roughly  correspond  to  the  vapor  pressure  at  100  K,  Pv~10-^ 
torr  [1.4].  At  300  K  this  pressure  corresponds  to  [ H2O ] =5X10°  cm  .  For 
a  baked  and  pumped  metal  system  a  typical  outgassing  rate  is  10-^® 
torr •liter/sec*cm  [1.5].  Assuming  ^0  is  a  major  constituent  of  the 

outgas  and  using  design  dimensions  of  the  test  cell  with  a  period  of 
about  30  hours  between  the  end  of  pumping  and  the  end  of  the  experiment, 
Pu  ^=10-10(torr*Ht/sec*cm2)  *60  cm2  •  lO-’sec/ 1  lit=6X10-l!t  torr.  At  300  K. 
this  pressure  corresponds  to  [H20]=2xlOl  cm  .  For  one  atmosphere  of 
this  H2O  concentration  corresponds  to  an  impurity  level  of  one  part  per 
million.  However,  the  water  present  will  be  very  effectively  scavenged 
by  reducing  reactions  with  Cs  vapor.  A  reduction  in  ^0  concentration  by 
a  factor  of  300  has  been  acheived  by  introducing  Cs  vapor  into  a  10~^ 
torr  vacuum  [1.6].  In  light  of  scavenging  [ H2O ] = 1 0 ^  ^  cm-^  appears  to  be 
a  conservative  worst  case  estimate. 

Even  assuming  an  amount  of  water  vapor  impurity  equal  to  this  worst 
case  estimate,  it  can  be  shown  that  ion-mole.cule  chemistry  of  water 
vapor  in  the  plasma  will  not  interfere  with  the  primary  experimental 
goal  of  measuring  the  rate  of  reaction  (1.2).  The  association  and 
switching  reactions  (1.3)  and  (1.4)  can  occur  with  replaced  by  H2O; 
related  bond  energies  are  ECs+*N2  -6.7  kcal/mol  [1.7],  and  ECg-H20=13.7 
kcal/mol  [1.8].  In  relation  to  the  measurement  of  the  rate  of  reaction 
(1.2),  reactions  (1.3)  and  (1.4)  provide  a  channel  for  the  production  of 
Cs?  (which  subsequently  dissociat Ively  recombines)  and  a  loss  mechanism 
for  Cs  +  .  It  has  been  shown  for  the  relatively  fast  reactions  involving 
the  Cs+,N2  cluster  that  a  partial  equilibrium  calculation  can  be  used  to 
determine  species  concentrations  before  recombination  begins.  Partial 


equilibrium  represents  a  worst  case  since  it  allows  the  greatest 
advancement  down  the  reaction  channel  which  produces  Cs^  prior  to  the 
onset  of  recombination.  Partial  equilibrium  may  therefore  be  conserva¬ 
tively  assumed  for  the  reactions  involving  the  Cs+»H90  cluster  in  a 
calculation  which  is  analogous  to  the  Cs"*"*^  partial  equilibrium  compu¬ 
tation.  Agreement  between  analytical  and  numerical  solutions  for  Cs^.H^O 
concentration  in  Fig.  1.4  supports  the  valididty  of  the  partial  equilib¬ 
rium  assumption.  For  equilibrium  involving  Cs  +  ^l^O  the  governing  equa¬ 
tions  are  (1.10),  (1.12),  (1.13),  (1.14), 


tCs+-  H2°]  Sb 

[Cs+][H20]  K_3b 


[h2o] 


13  -3 

constant  =  10  cm 


(1.15) 


(1.16) 


Equation  (1.13)  still  accurately  decribes  charge  conservation  since 
from  equations  (1.15)  and  (1.16)  and  Taoxe  1.1,  [Cs^l^O]  /  [Cs+]  = 
[H2O] *k3b/k-3b  ^  0.01%  at  all  temperatures  greater  than  500  K.  Both  the 
Cs+*N2  and  Cs^f^O  clusters  have  been  shown  to  occur  in  such  small 
concentrations  (relative  to  Cs+)  that  the  Cs^  equilibrium  concentration 
can  be  determined  independently.  In  Fig.  1.4  the  values  for  the  Cs^  and 
Cs+  concentrations  have  been  calculated  first  and  then  used  to  determine 

^  X  a. 

the  Cs  *N2  and  Cs  ‘l^O  concentrations.  The  Cs  cluster  will  not 

provide  a  significant  loss  mechanism  for  Cs+  and  the  presence  of  water 
vapor  will  not  affect  the  equilibrium  concentration  of  Cs^  calculated 
from  Equations  ( l . 10) ,  (1.12),  (1.13)  and  (1.14).  Since  the  Cs?  concen¬ 
tration  is  unchanged,  previous  statements  about  the  negligible  effect  of 
dissociative  recombination  of  Cs?  still  hold  in  the  presence  of  water 

vapor  impurities. 

Conclusions 

A  model  has  been  constructed  to  assess  the  role  of  several  ion- 

molecule  reactions  [(1.3),  (1.4),  (1.5)  and  (1.6)]  as  well  as  two  other 
reactions  [(1.7)  and  (1.8)]  which  have  not  been  addressed  in  previous 

analyses  of  the  experiments  being  planned  to  measure  the  electron 


recombination  rate  in  plasmas  under  conditions  of  relatively  low 
electron  number  density  and  high  ambient  pressure.  The  results  of  this 
study,  based  on  numerical  and  analytical  solutions  of  a  full  set  of 
kinetic  equations  with  rate  coefficients  estimated  from  information  in 
the  literature,  indicate  that  for  temperatures  in  the  range  of  1000  K  to 
1500  K,  electron  decay  rate  measurements  during  the  first  ten  milli¬ 
seconds  of  plasma  relaxation  can  be  identified  with  the  three-body 
process  e  +  Cs+  +  N£  ♦  Cs  +  ^2»  an<^  can  therefore  be  used  to  test  theo¬ 
retical  predictions  for  this  process.  Our  analysis  of  possible  difficul¬ 
ties  associated  with  the  presence  of  a  water  vapor  impurity  indicates 
that  in  this  temperature  range  the  effect  of  the  impurity  is  negligble. 

It  should  be  noted  that  considerable  uncertainty  exists  for  the 
values  of  many  of  the  rate  constants  used  in  this  study.  However  the 
study's  conclusions  are  based  on  robust  arguments  about  partial  equili¬ 
bria  and  the  ordering  of  time  constants  which  generally  demonstrate  that 
two  reaction  rates  or  two  characteristic  times  differ  by  many  orders  of 
magnitude.  The  strength  of  the  inequalities  involved  make  these  conclu¬ 
sions  insensitive  to  reasonable  uncertainties  in  rate  constants. 

Experimental  Work 

Description  of  Pulsed  Photoionization  Experiment 

In  the  past  year  substantial  progress  has  been  made  toward  com¬ 
pleting  the  construction  of  the  facility  for  the  pulsed  photoionization 
experiment.  The  cesium  cell  with  ancillary  heating,  temperature  measure¬ 
ment,  and  vacuum  systems  has  been  completely  fabricated.  A  breadboard 
circuit  for  the  flashlamp  has  been  built  and  preliminary  pulse  shapes 
have  been  obtained.  This  section  describes  the  design  of  the  experiment. 

The  radiation  source  chosen  for  photoionization  is  a  pulsed  short 
arc  Xe  lamp  directed  along  the  axis  of  a  tube  containing  the  test  gas 
(e.g.  Cs  vapor  and  M9) .  A  schematic  of  the  arrangement  is  shown  in  Fig. 
1.7.  The  Cs  vapor  and  are  contained  in  a  one  inch  diameter,  twelve 
inch  long  inconel  tube  fitted  with  sapphire  windows.  The  Cs  pressure  is 
one  to  ten  torr  and  the  nitrogen  pressure  is  about  one  atmosphere.  An 


electric  tube  oven  heats  the  Cs-^  mixture  in  the  temperature  range  400 
to  1300  K.  The  sapphire  window  has  75%  transmission  of  radiation  down  to 
200  nm  wavelengths  so  it  is  transparent  to  photons  capable  of  singly 
ionizing  ground  state  Cs  atoms,  i.e.  X  <  318  nra,  and  opaque  to  photons 
capable  of  ionizing  N2,  i.e.  X  <  200  nm.  The  Xe  flashlamp  provides  a  2 
ms  pulse  of  ionizing  photons  which  is  coupled  to  the  Cs  test  cell  by  an 
ellipsoidal  reflector. 

A  nonequilibrium,  isothermal  Cs  plasma  is  produced  by  the  radiation 
pulse.  The  free  electrons  are  rapidly  thermalized  through  collisions 
with  the  N2.  As  the  electrons  begin  to  recombine  with  Cs+  in  the 
presence  of  the  nitrogen,  the  electron  number  density  will  be  determined 
as  a  function  of  time.  Methods  for  obtaining  the  electron  density 
include  measurements  of  the  spectral  emission  from  free-bound  and  high 
lying  bound-bound  transitions  in  Cs ,  as  well  as  measurements  of  the 
electrical  conductivity  using  a  four  pin  conductivity  probe.  The  time 
history  of  the  electron  number  density  is  directly  related  to  the  rate 
coefficient  for  recombination  on 

The  Cesium  Cell 

Inconel  601  (Huntington  Alloys)  was  chosen  as  the  Cs  tube  wall 
material  because  of  its  exceptional  ability  to  resist  oxidation  and 
spalling  at  temperatures  up  to  1530  K.  The  1  inch  diameter  tube  will  be 
heated  at  atmospheric  pressure  in  a  tube  oven  to  temperatures  in  the 
range  between  400  and  1300  K,  and  it  must  withstand  oxidation  by  air 
externally,  and  reduction  by  Cs  internally.  Where  the  inconel  tube 
extends  beyond  the  oven,  the  temperature  will  be  maintained  at  700  K  by 
two  independent  ovens.  This  end  temperature  is  the  maximum  working 
temperature  of  the  sapphire  windows  (necessary  for  transmitting  ionizing 
and  diagnostic  radiation),  instrumentation  feed-through  (for  the  conduc¬ 
tivity  probe),  and  vacuum  flanges.  Because  this  hardware  cannot  sustain 
the  entire  temperature  range  of  the  experiment,  the  cesium  cell  has  end 
regions  held  at  a  lower  temperature  than  the  center  region. 

The  Cs  vapor  pressure  will  be  controlled  by  carefully  monitoring  the 
temperature  of  the  Cs  reservoir.  The  reservoir  will  be  maintained  20  to 


50  K  cooler  chan  any  other  part  of  the  system  to  prevent  Cs  condensation 
In  other  parts  of  the  system.  The  oxygen-free  copper  reservoir  and  the 
N2  supply  and  vacuum  pump  are  attached  to  the  Cs  cell  by  stainless  steel 
tubing.  The  diffusion  of  the  Cs  vapor  from  the  reservoir  through  the  one 
inch  diameter  by  12  inch  long  tube  to  90%  of  final  concentration  will 
occur  in  the  time  t*0.47  x^/D  [1.24].  Here  x  is  the  tube  length  and 
D=0.0911  (T/  288)^2  cm^/sec  is  the  ^-Cs  binary  diffusion  coefficient. 
For  T=*400  K,  t  =  1 .85  hours.  Although  not  unacceptably  long,  this  mixing 
time  will  be  decreased  by  natural  convection.  The  option  exists  to 
spread  the  Cs  throughout  the  reaction  cell  by  flowing  N2  through  the 
cell  at  a  low  flow  rate  relative  to  the  rate  of  chemical  reactions.  The 
flowing  system  also  has  the  advantage  of  constantly  removing  outgassed 
impurities  from  the  reaction  cell. 

The  1^2  number  density  will  be  determined  from  simultaneous  temper¬ 
ature  and  pressure  measurements  before  the  cell  is  heated.  To  not  exceed 
the  1%  uncertainty  estimated  for  the  thermocouple  measurement,  a  1% 
limit  on  pressure  measurement  uncertainty  is  set.  This  accuracy  must 
hold  over  the  pressure  range  200  to  1500  torr  (resulting  from  a  final 
pressure  of  1500  torr  with  final  temperatures  ranging  from  400  to 
1300K.)  These  requirements  will  be  met  by  a  series  of  Wallace  and 
Tiernan  gages  which  span  the  desired  pressure  range. 

The  tube  and  reservoir  is  attached  via  a  bakeable  valve  to  a  liquid 
nitrogen  trapped  diffusion  pump.  Pressures  of  lO-^1  torr  are  achievable 
with  this  system  in  order  to  remove  impurities  during  initial  bake  out 
and  subsequent  gas  changes. 

Photolonlzatlon 

Photoionization  of  the  Cs  by  Xe  flashlamp  radiation  occurs  through 
both  one  and  two  photon  processes.  The  single  photon  process  converts  a 
ground  state  6s  electron  to  a  free  electron  as  represented  by  the  reac- 


2  -f 

Cs(6  +  +  Cs  +  e 


26 


Here  the  photon  energy,  hv^,  must  exceed  the  ionization  potential  of  Cs , 
i.e.  X< 3 1 8  nm.  Two-photon  ionization  first  excites  a  ground  state  elec¬ 
tron  with  a  resonant  photon;  a  second  photon  then  ionizes  the  excited 
atom.  The  two  steps  are 

Cs(62S1y2)  +  hv2  -*•  Cs  (62P ^2) 

Cs*(62p3/2>i/2)  +  hv3  >  Cs+  +  e 

The  resonance  wavelengths  are  *2=852  nm  (J=3/2)  or  *2=894  nm  (J=l/2) , 

and  ionization  from  the  6p  state  can  occur  with  any  photon  of  wavelength 

shorter  than  *3<508  nm  (J=3/2)  or  *3<494  nm  (J=l/2).  The  cross  sections 
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for  ionization  from  the  6s  and  6p  levels  are  approximately  10  cm  and 
10  cm  respectively.  Contributions  from  the  two  photon  process  to 
total  free  electron  production  are  estimated  to  be  substantial  [1.25]. 

The  experimental  arrangement  is  greatly  simplified  by  positioning 
the  light  source  away  from  the  Cs  cell  and  the  heat  of  the  reaction 
zone.  Placing  the  source  outside  the  cell  necessitates  directing  the 
light  through  a  sapphire  window  at  the  end  of  the  inconel  tube.  To 
produce  a  collimated,  intense  photon  flux,  a  high  radiance  lamp  is 
required.  We  have  chosen  a  five  atmosphere  Xenon  short  arc  lamp,  (L- 
2369,  manufactured  by  ILC  Technology  of  Sunnyvale,  CA.)  This  lamp  has  a 
seven  mm  arc  length  and  a  synthetic  fused  silica  envelope.  It  is 
designed  for  pulsed  operation  at  a  maximum  energy  input  of  300  Joules 
per  pulse.  Pulsed  operation  [1.26,  1.27]  and  the  synthetic  fused  silica 
envelope  enhance  ultraviolet  output  from  the  lamp.  The  short  arc  lamp 
will  be  aligned  along  the  axis  of  an  ellipsoidal,  aluminum  reflector. 
(The  A-5000,  1000  watt  reflector  is  sold  by  Photon  Technology  inc., 

Princeton,  N.J.).  The  arc  is  imaged  at  the  far  focus  with  f/2.5  optics 
by  positioning  the  arc  at  the  near  focus  of  the  ellipsoid.  The  resulting 
beam  diverges  with  a  half  angle  of  11  degrees,  and  it  can  be  aimed 
directly  through  the  sapphire  window  into  the  tube.  Grazing  incidence 
and  high  reflectivity  of  the  polished  metal  tube  interior  surface  will 
insure  a  uniform  photon  flux  through  the  cell's  length. 


A  multiple  mesh  pulse  forming  network  has  been  prototyped;  it  pro¬ 
vides  a  roughly  rectangular  current  pulse  to  the  arc  lamp.  Pulse  width 
is  3  ms,  and  pulse  height  will  be  2,000  amps.  To  insure  negligible 
interference  by  photoionization  during  the  measured  plasma  decay,  the 
current  pulse  falls  from  its  maximum  value  to  zero  in  0.9  ms. 


Diagnostics 

A  four  pin  plasma  conductivity  probe,  developed  at  Stanford  in 
connection  with  research  on  MHD  plasmas,  appears  to  be  well  suited 
for  the  work  proposed  here.  The  principle  of  operation  of  the  probe 
is  to  measure  the  current  I  that  flows  when  an  ac  voltage  is  applied 
between  the  outer  pair  of  pins  (i.e.  electrodes)  of  a  linear  array  of 
four  pins,  and  to  measure  simultaneously  the  resulting  potential 
difference  6V  between  the  inner  pair  of  pins.  Since  J  (current  density) 
is  proportional  to  I  and  <5V  is  proportional  to  E  (electric  field),  it 
follows  from  Ohm's  Law,  E  =  J/o,  that  <5V  -  K  I/u  and  therefore  that 
o  »  K  dI/d(5V).  Here  o  denotes  the  electrical  conductivity  and  K  is  a 
calibration  constant  that  depends  only  on  the  geometric  features  of  the 
probe.  (In  practice,  5V  »  K  I/o  +  constant,  but  this  constant  plays  no 
role  upon  differentiation.)  The  constant  K  can  be  determined  by 
immersing  the  probe  in  a  liquid  electrolyte  (e.g.  aqueous  K.C1)  solution. 

The  conductivity  is  then  related  to  the  electron  number  density  by 
the  relation  a  =•  (nge  )/(meve^) ,  where  is  the  electron-heavy  colli¬ 
sion  frequency  and  e  is  the  electronic  charge.  Temporal  resolution  of 
50ys,  corresponding  to  operation  at  20  kHz  has  been  demonstrated 
[1.28].  Initially  this  probe  will  be  tested  on  a  steady  state  thermal 
plasma.  An  equilibrium  electron  concentration  of  3.6  x  10^  cm-^  can  be 
achieved  by  heating  the  Cs  cell  to  1300  K.  and  the  Cs  reservoir  to  630  K 
corresponding  to  a  vapor  pressure  of  5  torr.  The  four  pin  probe  measure¬ 
ment  will  be  compared  to  equilibrium  calculations  of  electron  concentra¬ 
tion  based  on  Cs  reservoir  and  cell  temperature  measurements.  The  probe 
will  also  be  compared  to  optical  diagnostic  measurements. 

Measurements  of  the  recombination  radiation  continuum  have  been  used 
to  measure  electron  densities  as  low  as  3X10^  cm"^  in  the  afterglow  of 


a  Cs-Ar  discharge  [1.29,  1.30].  In  this  work.  [1.29,  1.30]  a  monochro¬ 
mator  was  used  for  spectral  filtering.  By  employing  a  50  A  filter  in 
place  of  the  monochromator,  calculations  indicate  it  should  be  possible 

to  increase  the  collected  radiation  by  a  factor  of  the  order  of  10^  and 
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therefore  to  lower  the  limit  of  measureable,  transient  ng  to  about  10 
cm'\  This  limit  can  be  lowered  still  further  to  10^  cra"^  when  measure¬ 
ments  are  made  on  a  steady  state  plasma.  Hence  optical  emission  diag¬ 
nostics  should  provide  a  check  for  the  conductivity  probe  during  testing 
on  the  thermal  plasma.  Upon  successful!  testing  in  a  steady  state 
plasma,  the  conductivity  probe  will  be  used  to  measure  electron  densi¬ 
ties  in  the  plasma  afterglow. 

At  present,  a  simulated  probe  and  high  impedance  amplifier  have  been 
built  and  tested  with  a  1  MO  simulated  plasma  resistance.  Successful 
performance  has  been  obtained  at  frequencies  up  to  500  Hz.  Higher  fre¬ 
quency  operation  is  hampered  by  capacitive  coupling  between  probe  leads. 
The  problem  of  stray  capacitance  is  currently  being  considered. 


Sunmai 


A  detailed  analytical  and  numerical  model  has  been  constructed  and 
used  to  gain  understanding  of  reaction  kinetics  in  the  pulsed  photoioni¬ 
zation  experiment.  The  model  has  demonstrated  the  feasibiltiy  of  deter¬ 
mining  the  three-body  recombination  rate  of  the  reaction  e  +  Cs+  +  N2  + 
Cs  +  ^2  through  measurement  of  the  rate  of  electron  concentration  decay, 
provided  expriments  are  performed  at  temperatures  of  approximately  1000K 
or  greater.  Much  of  the  experimental  apparatus  has  been  constructed.  The 
reaction  cell  and  ancillary  heating  and  vacuum  systems  have  been  com¬ 
pleted,  and  extensive  development  has  been  done  on  the  pulsed  flashlamp 
system. 


2.2  Discharge  Effects;  Plasma-Electrode  Interaction 


Introduction 

The  structure  of  the  plasma-sheath  electrical  boundary  layer  between 
an  electrode  and  a  weakly  ionized  thermal  plasma,  and  the  mechanisms  of 
current  transfer  are  subjects  of  considerable  importance  for  a  number  of 
plasma  devices  and  applications.  These  include  arc  discharges, 
thermionic  converters,  MHD  generators  and  MPD  thrusters.  Important 
practical  aspects  or  the  interaction  of  electrodes  with  a  plasma  concern 
the  voltage  drop  between  the  plasma  and  the  electrode,  and  the 
conditions  determining  the  transition  from  diffuse  to  constricted  modes 
of  current  transfer.  Operation  of  both  cathodes  and  anodes  in  a 
constricted  mode  leads  to  electrode  erosion  which  sets  a  practical  limit 
to  electrode  life. 

While  it  has  been  the  subject  of  considerable  research  over  many 
decades,  the  plasma-electrode  boundary  layer  is  still  poorly  understood 
[2. 1,2.2].  The  difficulty  in  analyzing  this  region  stems  from  the 
inherent  nonlinearities  which  result  from  the  interaction  of  electrical, 
thermal,  fluid  dynamic  and  electrode  surface  effects. 

This  work  is  directed  towards  providing  a  more  complete 
understanding  of  the  plasma-electrode  interaction  for  the  case  of 
flowing,  seeded,  thermal  plasmas  at  pressures  of  the  order  of  one  atmo¬ 
sphere,  and  involves  both  theoretical  and  experimental  components. 

Consider  a  channel  flow  of  thermal  plasma  which  might  be  either  a 
noble  gas  or  combustion  products  gas  seeded  to  a  suitable  level  (~  1%) 
with  cesium  or  potassium.  We  may  define  the  electric  boundary  layer  as 
that  region  near  the  channel  walls,  composed  of  high  temperature  insula¬ 
tors  and  electrodes,  where  charged  particle  nonequilibrium  effects 
become  significant  and  even  dominant.  Normally,  the  walls  have  to  be 
operated  at  temperatures  less  than  that  of  the  core  flow  plasma  because 
of  material  limitations. 

On  physical  and  dimensional  grounds  the  region  between  the  core 
plasma  and  the  electrode  surface  can  be  divided  into  a  number  of  over¬ 
lapping  layers  of  decreasing  scale.  Proceeding  from  the  core  plasma  to 
the  electrode  these  layers  may  be  listed  as  follows  (see  Fig.  2.1): 


The  gasdynamic  boundary  layer  of  scale  5  ~  10  m, 
depending  on  distance  along  the  channel,  and  com¬ 
prising  both  a  velocity  and  thermal  boundary  layer. 
In  this  region  the  plasma  is  in  thermodynamic 
equilibrium  at  the  local  temperature  T(x)  and  supports 
an  electric  field  E(x)  =  J/a(x)  where  J  is  the  current 
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(a)  J»0 


Figure  2.1 


Potential  distribution  across  plasma  diode. 

I  Equilibrium  Plasma 

II  Ionization  Non-Equilibrium  Layer 

III  Sheath.  (Not  to  scale) 


density  and  a(T)  is  the  equilibrium  electrical  con¬ 
ductivity. 

)  The  ionization  nonequilibrium  (or  ambipolar  diffusion) 

layer  of  scale  TR  =  (2Da/Bn£)^^  ~  10~^m.  Here 

A  (T)  ~  10-2®  T-^/ 2  m*Vsec  is  the  three-body  recombina¬ 
tion  coefficient  divided  by  the  plasma  concentration 
at  infinity,  n^,  and  Da  is  the  ambipolar  diffusion 
coefficient.  In  this  layer  the  generation  rate  of 
electron-ion  pairs  exceeds  the  recombination  rate  to 
supply  net  fluxes  of  charged  particles.  In  the  absence 
of  current  (J=0)  the  equal  electron  and  ion  fluxes 
flow  to  the  electrode  where  they  recombine.  When  J*0, 
and  in  the  absence  of  electrode  emission,  these  fluxes 
flow  towards  or  away  from  the  electrode  and  constitute 
the  only  source  of  net  current.  A  potential  drop  of 
scale  (kT/e)  develops  such  that  the  walls  are  normally 
negative  relative  to  the  plasma,  except  for  an  anode 
at  large  currents  when  the  potential  fall  becomes 
positive  (see  broken  line  in  Fig.  2.1).  The  self- 


consistent  field  adjusts  the  charged  particle  fluxes 
to  match  the  core  current  to  that  entering  the  sheath 
overlying  the  electrode.  In  this  layer  quasi- 


neutrality  n  '  n.,  exists. 


(iii)  The  sheath  of  scale  equal  to  the  Debye  length  XQ  = 
(e^T/n^e  typically  lCf^ra.  In  this  layer  charge 

neutrality  is  not  preserved  and  large  electric  fields 
develop  which  must  be  calculated  from  Poisson's  equa¬ 
tion.  The  sheath  serves  to  match  the  charged  particle 
fluxes  from  the  ambipolar  diffusion  layer  to  the 
absorption  and  emission  processes  for  charged  parti¬ 
cles  at  the  electrode  surface.  For  a  cathode,  the 
sheath  structure  and  potential  drop  is  critically 
dependent  on  the  electron  emission  process  and  thus 
depends  on  the  temperature  and  surface  state  which 
control  the  work  function.  With  electron  emission,  a 
double  sheath  forms  at  a  cathode,  and  consists  of  an 
electron  rich  sheath  next  to  the  electrode  and  the 
usual  positive  ion  sheath  between  the  electron  sheath 
and  the  neutral  plasma.  For  an  anode  that  does  not 
emit  tons,  the  current  transfer  is  effected  purely  by 
the  collection  of  electrons  from  the  plasma.  The 
situation  is  different  in  those  cases  (e.g.,  in 
thermionic  converters)  where  the  anode  can  emit  ions. 


(iv)  The  Knudsen  layer  where  free  molecule  flow  takes 
over  from  a  continuum  (fluid)  description  and  the 
appropriate  scale  is  the  mean  free  path.  The  mean 

free  paths  for  electrons  and  ions  are  of  order 

4e  —  ~  10~^m  and  are  thus  comparable  to  Xp.  Our 

calculations  show  that  the  sheath  is  typically  many 
Debye  lengths  thick,  so  that  the  sheath  is,  for  the 
most  part,  collision  dominated. 


It  should  be  emphasized  that  the  four  layers  distinguished  above  are 
overlapping  and  no  sharp  boundaries  can  be  distinguished.  Moreover,  the 
scale  lengths  are  nominal  ones  based  on  core  values  rather  than  local 
conditions.  In  general  the  scale  lengths  are  ordered  according  to  5  >> 

eR  >>  V  D  ~  *e  ~  lf 

As  the  current  density  J  is  increased  from  zero,  the  potential  drop 

at  the  cathode  increases  and  more  ions  than  electrons  are  collected. 

2 

However  this  condition  only  holds  for  very  small  currents  (£  1  mA/ cm  ) , 
at  which  point  the  collected  electron  flux  becomes  zero  and  the  current 
at  the  cathode  is  carried  entirely  by  arriving  ions.  For  higher 


currents,  of  technical  interest,  the  cathode  must  emit  electrons.  Of  the 
known  emission  processes,  only  thermionic  emission  (including  field- 
assisted,  Schottky  emission)  is  relevant.  Depending  on  conditions, 
particularly  the  cathode  material,  temperature  and  work  function,  dif¬ 
fuse  emission  persists  up  to  some  current  density  when  there  is  a  tran¬ 
sition  to  a  constricted  discharge  terminating  in  one  or  more  arc  cathode 
spots . 

At  the  anode  it  might  be  expected  that  large  currents  up  to 
J  ~  en^  ce/4,  the  random  electron  current  (~  250  A/cm^),  could  be  col¬ 
lected  as  a  result  of  the  reducing  sheath  potential  drop.  However  it  is 
commonly  observed  that  constrictions  in  the  form  of  anode  glow  spots 
occur  at  much  lower  currents  and  that  the  sheath  potential  drop  changes 
sign  as  indicated  in  Fig.  2.1.  The  reason  for  this  behavior  is  hitherto 
unexplained,  but  it  is  generally  thought  to  be  related  to  the  creation 
of  a  highly  conducting  discharge  channel  which  penetrates  the  lower 
conductivity  boundary  layer  plasma. 

From  the  foregoing  discussion  it  is  clear  that  the  physical  pro¬ 
cesses  involved  are  extremely  complex  and  that  experimental  data  are 
vital  to  guide  the  theoretical  modeling. 

Research  Objectlvea 

t.  To  observe  and  record  the  current-voltage  characteristics  and  the 
transition  from  diffuse  to  constricted  discharge  modes  for  elec¬ 
trodes  immersed  in  a  flowing,  atmospheric  pressure  thermal  plasma. 

2.  By  analytical  and  computer  modeling,  to  predict  the  current-voltage 
characteristics  for  anodes  and  cathodes  under  diffuse  conditions 
and  to  compare  these  predictions  with  the  experimental  results. 


3. 


To  use  the  insight  provided  from  objectives  1  and  2  above  to  more 
fully  understand  the  constriction  process. 


Status  of  Research 


To  date,  we  have  concentrated  our  efforts  on  the  theoretical  aspects 
of  the  research  objectives. 

Theoretical  Work 

The  approach  adopted  for  theoretical  modeling  is  to  start  off  with 
the  simplest  relevant  model  and  to  introduce  the  complexities  of  the 
complete  problem  one  at  a  time.  To  this  end,  the  plasma-sheath  problem 
for  a  plane  electrode  in  contact  with  a  uniform,  equilibrium,  weakly  - 
ionized  plasma  has  been  formulated  in  some  generality.  By  assuming  that 
the  plasma  is  isothermal  and  stationary,  the  complexities  associated 
with  the  gasdynainic  thermal  boundary  layer  and  Joule  heating  are,  for 
the  present,  eliminated  to  allow  the  essential  plasma  physics  aspects  of 
the  problem  to  be  studied. 

This  formulation  in  terms  of  the  conservation  equations  foe  the 
electron  and  ion  concentrations  and  momenta  (including  inertia), 
together  with  Poisson's  equation,  gives  a  uniform  description  of  the 
ionization  nonequilibrium  layer  and  sheath,  in  both  the  collionless  and 
collision  dominated  regimes,  without  having  to  resort  to  an  arbitrary 
division  and  matching  of  the  separate  regions. 

To  the  present  time,  under  this  grant,  this  plasma-sheath  problem 
has  been  investigated  in  the  continuum  (collision-dominated)  limit, 
neglecting  the  inertia  terras,  which  is  applicable  to  MHD  and  similar 
high  pressure  plasma.  In  this  case  the  governing  equations  are: 


n  =  f  1  -  n  n .)  n 
v  e  v  i 
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together  with  the  subsidiary  relations 


(2.5) 


J  =  - in —  j-r  _  r  j 

(  1  +  U  1  6 


(2.6) 


<j>  =  -  f  E  dx 


(2.7) 


?  1/2 

Here  x  is  normalized  on  the  ion  recombination  distance  In  =  (2D  /Bn  )  ' 

k  a  °° 

where  Da  is  the  ambipolar  diffusion  coefficient  and  S  is  the  three-body 
recombination  coefficient.  The  electron  and  ion  concentrations  are 
normalized  on  n  while  the  net  generation  rate  n  is  normalized  on  3  nJ 

OO  CO  • 

The  potential  4>  is  normalized  on  (kT/e)  while  the  electric  field  is 
normalized  on  Er  =  (kT/e5,R).  The  fluxes  are  normalized  on  = 

(2D  n  /i  )  and  the  current  density  is  normalized  on  J  =  en  (u  +  u  •  ) 
En  =  o  E„  where  a  is  the  electrical  conductivity  for  x  +  <=». 

IV  OO  ^  30 

There  are  three  parameters  to  be  specified:  the  mobility  ratio 
U  =  (iij7ua)  «  1,  a  function  only  of  the  type  of  ion  and  background  gas; 
S  =  '^D^R^  <<  »  t-^le  rati°  °£  the  Debye  length  to  the  ion  recombination 

length;  and  J,  the  normalized  current  density. 

For  a  perfectly  absorbing  (catalytic),  non-emitting  electrode,  tne 
boundary  conditions  at  x=0  are  the  familiar  continuum  ones: 


n  (U)  =  n  (0)  =0 
e  i 


(2.8) 


For  an  electrode  which  emits  electrons  (or  ions)  the  continuum  boundary 

condition  is  that  n  (0)  (or  n.(0))  has  a  prescribed  non-zero  value. 

e  l 

However,  it  is  necessary  to  relate  the  value  of  ng(0)  (or  (n^(0))  to  the 
actual  value  of  the  emission  flux  density  from  the  surface.  By  con¬ 
sidering  the  evolution  of  the  charged  particle  distribution  function 
across  the  Knudsen  layer,  the  following  relation  has  been  derived 


relating  the  electron  saturation  emissive  flux  J*(U)  at  the  surface  and 
the  values  of  n^CO)  and  ^(O)  in  the  continuum  description. 

T(0)  =  -  -2jJ--2-  [n  (0)  (ce/2v  )  +  r  (0)]  (2.9) 

(1  +  V) 

where  vD  =  (2D  /$.  )  is  a  reference  diffusion  velocity  and  c  = 

&  .  a  R  e 

1/2 

(8kT/vme)  ”  is  the  elec  tron  mean  thermal  speed. 

Physically,  J*(0)  is  the  electron  emission  saturation  current 
density  given  by  the  Richardson-Dushman  law 

Je(0)  =  -  AT  exp  -  Ow/kT)  (2.10) 

c  2  2 

where  <J>W  is  the  surface  work,  function  and  A  =  1.2  x  10°  A/m  K  .  Equation 
(2.10)  is  applicable  when  the  electric  field  at  the  electrode  is  zero  or 
decelerating  for  emitted  electrons  (E(0)  >  0).  When  there  is  an 

accelerating  field  at  the  electrode  (E(0)  <  0),  field-enhanced  emission 
due  to  the  Schottky  effect  is  applicable; 

J>)  =  -  AT2  exp  -  [$w  -  (-  e3E(0)/4ireo)  ^2]/kT  (2.11) 

“  JgG  exp  (0.44  (-  E(0) )1/2/t)  (2.1lb) 

where  JgQ  is  the  emission  given  by  Eq.  (2.10).  In  practice  very  large 
accelerating  fields  are  required  to  cause  significant  enhancement  of  the 
emission . 

The  case  of  ion  emission  at  the  electrode  likewise  corresponds  in 
general  to  the  boundary  condition  n^(0)  >  0.  However,  ion  emission  only 
occurs  for  the  relatively  rare  combination  of  a  gas  whose  ionization 
energy  is  comparable  with  or  less  than  the  work  function  of  the  (bare) 
surface,  e.g.,  the  heavier  alkali  metal  vapors  combined  with  refractory 
metals,  such  as  are  employed  in  thermionic  converters.  Depending  on  the 
electrode  temperature,  the  surface  is  more  or  less  covered  with  a  layer 
of  the  alkali  metal,  which  lowers  the  work  function  and  generally  pro¬ 
motes  electron  emission  as  well  as  ion  emission.  For  most  plasma 


devices,  including  MHD  and  MPD  applications,  the  possibility  of  signifi¬ 
cant  ion  emission  is  negligible. 

The  remaining  boundary  conditions  for  the  equation  set  (2.1)-(2.5) 
are  taken  as  that  the  values  of  r^,  and  E  at  x  «  should  match  those 
given  by  the  quasi-neutral  approximation  [2.3]  which  will  now  be  briefly 
summarized . 

In  the  quasi-neutral  solution  we  put  n  =  n.  =  n  and  discard 

e_  1 

Poisson's  equation.  With  the  boundary  conditions  n(0)  =  0,  n(“)  =  1 ,  we 
obtain  the  analytical  results 


n  =  tanh  x 


,2  — 

n  =  sech  x  tanh  x 


T  =  -  7  sech2  x  -  J 

e  2  4p 


vi  ’  ~  2  sech2  x  +•  ^  J 


E  -  [-  ( |  +  u )  sech2  x  +  j]  coth  x 


(2.12) 

(2.13) 

(2.14) 

(2.15) 

(2.  In) 


-)>  »  £n  tanh  x  -  J  in  sinh  x 


It  is  noteworthy  n  and  n  are  independent  of  J,  and  that  E  has  a 
simple  pole  while  i  has  a  logarithmic  pole  at  x  =  U.  The  profiles 
of  n  and  n  are  shown  in  Fig.  2.2a,  while  the  potential  profiles  for 
various  J  are  shown  in  Fig.  2.2b. 

Within  the  past  year,  a  numerical  code  for  solution  of  the  lull 
equations  has  been  perfected.  This  uses  the  method  of  quas i 1 iaear izat ion 
and  orthonormalization  of  the  governing  equations.  A  large  number  of 
numerical  solutions  have  been  obtained  over  wide  ranges  of  the  parameter 
values  u,  e  and  J  with  different  levels  of  electron  and  ion 
emission,  corresponding  to  different  values  of  the  boundary  conditions 
ng(0),  n^(0).  Our  aim  is  to  provide  a  comprehensive  physical  under¬ 
standing  of  the  nature  of  the  various  forms  of  solution  and  to  document 


Figure  2.2a.  Profiles  of  normalized  charged 
particle  concentrations  and  net  generation 
rate . 
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Figure  2.2b.  Profiles  of  the  norraalized_potential  as 
function  of  normalized  current  density:  J  <  0,  cathod 
J  >  0,  anode.  Broken  lines  indicate  the  asymptotes 
corresponding  to  the  uniform  electric  field  in  the 
distant  plasma  for  the  quasi-neutral  approximation. 


the  dependence  of  various  important  quantities  such  as  the  sheath  thick¬ 
ness  and  potential  drop  as  a  function  of  the  parameters. 

The  results  have  been  organized  into  three  basic  cases:  (i)  an 

electrically  isolated  (floating)  electrode  (J=0)  without  and  with  elec¬ 
tron  emission,  (ii)  a  cathode  (J<0)  without  and  with  electron  emission 
and  (iii)  an  anode  (J> 0)  without  and  with  both  ion  and  electron  emis¬ 
sion.  The  first  case,  which  has  been  completed  and  written  up  [2.4]  is 

relevant  to  the  electrodes  of  MHD  generators  at  open  circuit  and  is  also 
related  to  the  theory  of  thermionically  emitting  probes,  first  intro¬ 
duced  by  Langmuir.  The  second  and  third  cases,  for  which  reports  are 
currently  being  written,  are  relevant  respectively  to  the  cathodes  of  a 
variety  of  high  pressure  plasma  devices  and  to  the  emitters  of  ther¬ 
mionic  converters  in  the  ignited  mode  as  well  as  to  the  anodes  of 

seeded,  noble  gas  MHD  and  MPD  devices.  A  few  of  the  important  results 

for  case  (i)  may  be  illustrated,  as  follows. 


for  J=0_ 

and 

zero  electron  emission  (ng(0)  - 

0), 

,  the  profiles 

of 

ne*  ni»  n 

,  and 

E  and  6 

are  illustrated  in  Figs. 

2.3 

and  2.4,  for 

the 

cases 

of  u 

=  1/300 

(typical  of  combustion 

MHD 

plasmas)  and 

e  =* 

o 

t 

o 

-2 

• 

In  Fig. 

2.3a,  for  c  =  10  *,  it 

is 

seen  that  both 

n  and  n,  lie  well  below  the  value  n_v,  given  by  the  quasi-neutral 
e  i  QN 

approximation,  with  n^  >  n^,  corresponding  to  a  positive  ion  sheath.  As 
a  result  the  position  of  maximum  net  generation  n  occurs  further 

from  the  electrode  (x  =»  1.125)  than  that  of  the  quasi-neutral  solution 
(x  =  0.658).  If  we  define  the  sheath  strength  as 

A  -  I  ni  "  ne  I  1  tni  +  neV2  (2.18) 

and  (arbitrarily)  define  the  plasma-sheath  boundary  as  corresponding 
to  A  =  0.1,  then  the  boundary  occurs  at  x  =  0.815,  corresponding  to 
a  sheath  thickness  of  8.15  Ag.  As  e  is  decreased  (Fig.  2.4a)  the  sheath 
is  thinner  on  the  x  scale,  but  on  the  Debye  length  scale  it  increases 
to  17.14  Aq*  Also,  as  e  is  decreased,  the  position  of  maximum  genera¬ 
tion  tends  to  that  given  by  the  quasi-neutral  solution.  The  sheath 
thickness  and  position  of  maximum  generation  are  given  as  a  function  of 


Figure  2.3  Profiles  of:  a)  n^,  ,  n  compared  with  those  of  quasi 

neutral  approximation,  n^,  n  ;  b)  E  and  parameter  values 
U  =  1/300,  e  =■  10~l,  7  =  0,  n  (0)  =  n  (0)  =  0. 
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e  (10  10  )  in  Table  2.1.  It  may  be  noted  that  for  typical 

_3 

combustion  MHD  conditions  (e  =  10  )  the  sheath  is  36.3  Debye  lengths 

thick.  Since,  for  the  same  conditions,  the  electron  and  ion  mean  free 
paths  are  comparable  to  the  Debye  length  (~  10  ^m),  this  justifies  the 
assumption  of  a  collision-dominated  sheath  inherent  to  the  continuum 
analysis.  The  profiles  of  E  and  <J>  in  Figs.  2.3b  and  2.4b  show 

that  E  increases  strongly  and  4>  increases  moderately  as  £  decreases.  The 
values  of  the  surface  field  E(0),  and  the  total  potential 

drop  are  given  in  Table  2.1  as  a  function  of  £.  Another 

interesting  result  is  that  7^(0)  =  I\(0),  which  equals  the  integral 
of  the  net  generation  rate  n,  remains  fixed  at  the  value  of  -  1/2  gi^ven 
by  the  quasi-neutral  solution,  despite  the  change  in  the  profile  of  n  in 
the  presence  of  the  sheath. 

Table  2.1  Non-emitting  Electrode  (J  =  0)  :  Ion  Sheath  Thickness  (in 
Debye  lengths),  Position  of  Maximum  Generation  Rate, 
Surface  Field,  and  Total  Potential  Drop  vs  £  for 
U  =  1/300. 


Position 

of 

£ 

Ax  /X 
is  D 

• 

n 

max 

E(0) 

A$ 

,  „  ./CkT/e) 
total 

10-1 

8.51 

1.13 

-13. 2 

-7.53 

10_3/Z 

11.73 

0.86 

-28.5 

-8.30 

10~2 

17.14 

0.76 

-bl  .3 

-9.06 

10-5/2 

25.20 

0.69 

-132.1 

-9.83 

to-3 

36.30 

0.68 

-284.5 

-10.61 

The 

profiles  of  n  , 

—  • 

n^,  n,  and  of 

E  and  j>  for  u 

1/300  and  e  =  10  1 

with  increasing  values 

of  electron 

emission,  ng(0) 

= 

0.1,  0.5  and  1.0 

are  shown  in  Figs.  2.5 

through 

2.7. 

The  presence 

of 

electron  emission 

Introduces  an  electron-rich  sheath  adjacent  to  the  electrode,  followed 
by  the  usual  ion  sheath.  With  increasing  electron  emission,  the  width 
and  strength  of  the  electron  sheath  increases  while  the  strength  and 
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width  of  the  ion  sheath  decreases.  As  a  result,  the  total  potential 

drop  decreases  dramatically;  the  electrode  potential  floats  up  towards 

the  plasma  potential,  becoming  equal  to  it  for  ng(0)  =  1.  This 

behavior  is  consistent  with  the  use  by  Langmuir  and  others  of  a 

thermionically  emitting  probe  to  measure  the  plasma  potential.  For  a 

given  value  of  e,  there  is  a  critical  value  of  electron  emission, 

measured  by  ne(0)crit»  eclual  t0  0.223  for  e  =  10~^,  at  which  the  surface 

field  changes  sign  from  negative  to  positive  and  a  potential  minimum 

occurs  at  the  electrode  surface.  As  n  (0)  is  increased  above  the 

e 

critical  value,  the  potential  minimum  moves  out  from  the  electrode 

surface  and  the  depth  of  the  potential  well  increases,  as  is  clear  from 

Figs.  2.6b  and  2.7b.  There  is  also  a  minimum  in  n  associated  with  the 

e 

potential  minimum,  due  to  the  fact  that  the  electrons  approximately 
satisfy  the  Boltzmann  relation  n  «  exp  It  should  also  be  noted  that 
although  the  electrode  potential  floats  up  to  the  plasma  potential 
for  n^CO)  =  1,  the  plasma-sheath  region  is  not,  by  any  means,  a  region 
of  uniform  potential  and  charged  particle  concentration. 

To  our  knowledge,  this  work  represents  the  first  satisfactory  quan¬ 
titative  analysis  of  an  emitting,  floating  electrode,  and  shows 

significant  differences  from  the  very  approximate  treatment  originally 
given  bv  Langmuir.  Figure  2.8  shows  the  total  potential  drop  as  a 

function  of  electron  emission  fn  (0)1  for  various  values  of  e.  It 

v  e  ' 

indicates  that  except  for  very  small  values  of  emission,  the  potential 
drop  is  independent  of  e,  and  that  while  the  electrode  potential 
approaches  the  plasma  potential  for  relatively  small  values  of  ng(0),  it 
only  tends  to  exact  equality  as  ng(0)  tends  to  unity. 

One  other,  unexpected,  result  should  be  mentioned.  With  increasing 

electron  emission,  there  is  a  significant  shift  of  the  region  of  net 

generation  to  the  electrode  surface,  because  of  the  enhanced  value 

of  n  in  the  electron  sheath.  However,  despite  the  dramatic  change  in 

•  - 

the  prof ile^o^  n,  the  integral  of  the  net  generation  rate  (for  J =0)  i.e. 

Tt(0)  =  -  /  n  dx,  remains  essentially  invariant  at  the  value  -  1/2 

given  by  the  quasi-neutral  approximation. 
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Figure  2.8.  Total  potential  drop  &$total  as  a  functi-on  °f  ne(0),  tor 
various  t;  parameter  values,  u  =  1/300,  J  =  0,  n^(0)  =  0. 


Experimental  Work 


The  experimental  work  in  progress  is  based  on  the  use  of  a  versatile 
bench-scale  diffusion  burner  capable  of  producing  a  combustion  plasma  to 
JOOGK.  It  is  possible  to  seed  the  combustion  products  with  1%  by  mass  of 
potassium  vapor  to  produce  a  thermal  plasma  with  a  conductivity  of  up  to 
10  S/m.  This  facility  was  built  in  connection  with  an  earlier  h'SF  grant 
for  the  testing  of  plasma  diagnostic  techniques. 

During  the  past  contract  year  we  have  installed  and  partially 
configured  a  dedicated  personal  computer  to  be  used  for  data 
acquisition.  This  machine  (an  Intel  System  310)  was  donated  by  the  Intel 
Corporation  as  part  of  an  academic  grant  program. 

We  expect  to  measure  the  current-voltage  characteristics  and  perform 
breakdown  experiments  during  the  current  contract  year. 


2.3  Interaction  of  Discharges  and  Fluid 


.amics 


Introduction 

In  several  space  power  and  propulsion  systems  of  potential  interest 
to  the  Air  Force,  such  as  MPD  thrusters,  MHD  generators,  and  rail  guns, 
fluid  motions  induced  by  discharges  are  of  critical  importance  to  both 
the  understanding  and  performance  of  the  system.  The  channel  flow  of  a 
plasma  through  a  transverse  magnetic  field  is  influenced  by  the  Hall 
effect,  which  causes  nonuniform  Lorentz  forces  to  arise  in  the  direction 
mutually  perpendicular  to  the  channel  axis  and  the  magnetic  field.  These 
forces  cause  secondary  flows  to  develop  as  the  flow  progresses  down  the 
region  of  electromagnetic-fluid  interaction.  The  theoretical  treatment 
of  this  problem  has  included  the  work  of  Fay  [3.7],  Broer,  Peletier  and 
Wijngaarden  [3.4],  Sato  [3.18],  Tani  [3.20],  Apollonskil  and  Kos'kin 
[3.1],  Sastry  and  Bhadram  [3.17],  Liu,  Lineberry  and  Schmidt  [3.11],  and 
Ishikawa  and  Umoto  [3.9].  These  secondary  flows  were  predicted  to  be  of 
substantial  magnitude  in  the  calculations  reported  by  Maxwell  et  al 
[3.13],  Bityurin,  Zatelepin  and  Lyubimov  [3.3]  and  Doss  and  Ahluwalia 
[3.6].  The  calculated  magnitude  of  these  flows  range  from  10  to  30%  of 
the  bulk  velocity,  depending  on  conditions,  and  the  predicted  conse¬ 
quences  include  the  development  of  significant  asymmetries  in  the  pro¬ 
files  of  mean  axial  velocity  and  of  temperature,  with  resultant  effects 
on  wall  heat  transfer  and  electrode  voltage  drops.  In  addition 
Deraetriades  et  al  [3.3]  predicted  that  these  secondary  flows  could  lead 
to  raagnetoaerothermal  instabilities,  which  could  include  boundary  layer 
separation  and/or  electric  field  breakdown  near  the  electrodes. 

Indirect  evidence  for  the  existence  of  MHD-induced  secondary  flow, 
reviewed  by  Girshick  and  Kruger  [3.8],  included  measurements  of  temper¬ 
ature  profiles  obtained  by  James  and  Kruger  [3.10]  and  of  electrode 
voltage  drops  obtained  by  Barton  [3.2],  which  indicated  that  MHD  effects 
caused  the  anode  wall  boundary  layer  of  the  channel  to  be  colder  than 
expected  while  the  cathode  wall  is  hotter  than  expected.  It  is  hypothe¬ 
sized  that  these  results  could  be  ascribed  to  a  secondary  flow  mecha¬ 
nism.  Denetriades  et  al  [3.5]  conjectured  that  evident  arc  damage  to  the 
center  of  the  anode  wall  in  the  downstream  half  of  the  channel  used  in 
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Che  High  Performance  Demonstration  Experiment  at  the  ArnoLd  Engineering 
Development  Center  (Starr  et  al,  3.19)  is  a  manifestation  of  the  raagne- 
toaerothermal  instability  mentioned  above.  McClaine,  Swallom  and  Kessler 
[3.14]  reported  axial  velocity  profiles  measured  using  a  traversing 
stagnation  pressure  probe  mounted  in  the  diffuser  of  the  Avco  Mk  VI 
'  generator.  The  profile  measured  between  the  top  and  bottom  walls  is 

found  to  be  strongly  skewed  by  MHD  forces,  presumably  as  a  consequence 
of  secondary  flow. 

Our  current  experiments  provide  the  first  direct,  quantitative 
j  measurements  of  the  secondary  flow  field  in  an  MHD  channel.  Addi¬ 

tionally,  the  effects  of  secondary  flow  on  the  plasma  momentum,  thermal 
and  electrical  behavior  a^e  being  studied  by  making  measurements  of  the 
profiles  of  mean  axial  velocity,  turbulence  intensity  (axial  and  trans- 
j  verse)  and  plasma  voltage,  and  by  measuring  electrode  surface  temper¬ 

atures. 

Research  Objectives 

|  The  over-riding  goal  of  our  work  in  this  area  is  to  obtain  experi¬ 

mental  information  which  will  allow  us  to  describe  the  three-dimensional 
flow  field  in  an  MHD  device.  Specifically,  the  major  objectives  which 

;  presently  guide  this  work  are  the  following: 

! 

1.  To  study  the  effect  of  electromagnetic  body  forces  on  the  secondary 

flow  field.  What  are  the  magnitudes  of  these  secondary  flows,  and 

how  is  their  vortex  structure  affected  by  applied  magnetic  fields 

|  and  discharges? 

2.  To  study  the  effect  of  MHD  interaction  on  mean  axial  velocity  pro¬ 
files.  Are  they  affected  by  turbulence  suppression  caused  by  the 
magnetic  field,  and  to  what  extent  are  they  distorted  by  secondary 
flows? 

3.  To  study  the  effect  of  MHD  interaction  on  turbulence.  What  is  the 

|  detaiLud  nature  of  magnetic  field  turbulence  suppression,  including 
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Che  development  of  anisotropic  turbulence?  To  what  extent  are  turbu¬ 
lence  intensity  profiles  distorted  by  secondary  flows? 


Status  of  Research 

Theoretical  Background 

The  fluid  momentum  equation  for  an  MHD  plasma  can  be  written 

Du  — 

p  —  =  -  7p  +  V • t  +  J  x  B  (3.1) 

where  p  is  the  mass  density,  p  the  thermodynamic  pressure  and  x  the 
viscous  stress  tensor.  With  B  =  Bz,  the  presence  of  a  non-zero  Hall 
current  causes  the  J  *  B  Lorentz  force  to  have  a  component  in  the  y- 
direction.  Taking  the  curl  of  (3.1)  to  obtain  the  fluid  vorticity  equa¬ 
tion,  we  obtain  for  the  axial  vorticity  the  result 

DQx  ~  3JX 

(3.2) 

where  we  have  assumed  Bz  =  constant. 

It  is  thus  evident  that  the  nonuniformity  of  the  Hall  current 

density  between  the  core  and  the  sidewalls  introduces  a  rotational 

component  into  the  Lorentz  force.  For  the  first-order  profile  of  Jx 
pictured  we  would  expect  the  two-cell  flow  structure  shown  on  the  right 
of  Fig.  3.1  to  arise.  The  presumed  effect  of  such  a  secondary  flow  would 
be  to  transport  hot,  fast  core  plasma  toward  the  bottom  wall  while 
sweeping  relatively  cold,  slow  sidewall  fluid  into  the  boundary  layer  of 
the  top  wall. 

As  the  axial  vorticity  and  the  distribution  of  Jx  are  coupled 
through  (3.2),  the  actual  structure  of  the  secondary  flow  field  could  be 
considerably  more  complicated  than  the  simple  pattern  suggested  in  Fig. 
3.1;  indeed,  Jx  might  not  increase  monotonically  from  the  walls  to  the 

core,  in  which  case  the  axial  vorticity  could  change  sign  in  some 

regions  of  the  cross-plane.  It  should  also  be  pointed  out  that  Faraday 


Figure  3.1  Secondary  flow  mechanism  for  a  case  of  positive  net  Hall 
current . 

MHD  generators  would  suppress  the  flow  of  net  Hall  current  by  using 
segmented  electrodes  separated  by  insulators.  However,  the  Hall  effect 
would  still  cause  the  local  value  of  Jx  to  be  in  general  non-zero;  only 
the  integrated  value  of  Jx  over  the  cross-plane  would  vanish,  and  then 
only  in  the  absence  of  axial  leakage  currents. 

Description  of  experiaents 

Flow  Train  and  Run  Conditions 


MHD  secondary  flow  is  being  investigated  in  a  series  of  experiments 
conducted  at  the  High  Temperature  Gasdynamics  Laboratory,  with  the 
M-2  fLow  train,  shown  in  Fig.  3.2.  This  consists  of  a  combustor  with  a 


y 

Figure  3.2  M-2  flow  train. 


nominal  2-MW  thermal  input  rating,  a  large  mixing  plenum,  a  nozzle 
through  which  the  combustion  products  are  accelerated  into  a  duct  having 
a  square  cross-section  measuring  5.1  cm  on  a  side,  a  diffuser  and  an 
exhaust  system  including  a  scrubber  and  a  stack..  The  square  duct  con¬ 
sists  of  a  run-in  section  and  the  active  channel  itself,  which  is 
enclosed  in  a  2.6-T  water-cooled,  copper-coil  magnet. 

The  plasma  consists  of  the  combustion  products  of  ethyl  alcohol 
(C2H5OH)  burned  in  pure  oxygen  at  a  stoichiometry  of  1.05  (fuel  rich). 
The  fuel  is  seeded  with  potassium  hydroxide  (KOH),  such  that  potassium 
constituted  2.1%  of  the  total  products  by  mass.  The  total  flow  rate, 
54.4  g/s,  at  a  calculated  nozzle  exit  temperature  of  2750  K,  corre- 

A 

sponded  to  a  Reynolds  number  (based  on  hydraulic  diameter)  of  1.3  *  10  ; 
the  calculated  nozzle  exit  velocity  is  171  m/s.  This  value  has  been  used 
throughout  to  characterize  the  bulk  velocity  U^,  although  it  should  be 
noted  that  the  actual  bulk  velocity  decreases  somewhat  in  the  downstream 
direction  as  the  plasma  cools  and  becomes  more  dense.  The  flow  thus 
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being  turbulent,  the  boundary  layers  at  the  velocity  measurement  plane, 
which  is  21.4  hydraulic  diameters  downstream  of  the  nozzle  exit,  are 
fully  developed  (or  at  least  nearly  so)  in  the  absence  of  MHD  inter¬ 
action. 


Electrical  Configuration 

An  important  goal  of  the  experiments  has  been  to  isolate  the  effect 
on  the  flow  of  the  Hall  current,  since  the  Faraday  current  itself  has 
important  effects  because  it  induces  an  axial  Lorentz  force  given  by 
JyB.  It  is  also  considered  desirable  to  establish  conditions  which  allow 
insofar  as  possible  for  a  straightforward  interpretation  of  the  measured 
secondary  flow  field.  Accordingly  the  channel  is  connected  in  the  "Hall 
configuration",  with  a  pair  of  electrodes  at  the  upstream  end  connected 
through  an  external  load  to  a  pair  of  electrodes  at  the  downstream  end, 
as  shown  in  Fig.  3.3.  A  current  of  9.3  amperes  is  driven  through  the 
plasma  by  means  of  a  bank  of  batteries,  and  controlled  by  varying  the 
load  resistance.  Thus  the  experimental  situation  corresponds  to  the  case 
illustrated  in  Fig.  3.1,  with  a  positive  net  Hall  current. 

A  non-dimensional  parameter  which  indicates  the  strength  of  the 
electromagnetic-fluid  interaction  is  the  Stuart  number,  or  magnetic 
interaction  parameter,  S.  For  this  configuration  an  appropriate  defini¬ 
tion  of  S  is 


mU 

where  Ix  is  the  applied  Hall  current,  L  the  distance  from  the  first 
loaded  electrode  to  the  measurement  plane,  m  the  mass  flow  rate  and  U 
the  mean  axial  velocity.  The  interaction  parameter  in  this  form  repre¬ 
sents  the  ratio  of  the  Lorentz  twisting  force  to  the  axial  inertia  of 
the  flow.  The  value  of  S  in  the  present  experiments  is  about  l.l  at  the 
plane  of  the  velocity  measurements  and  slightly  less  at  the  measurement 
planes  for  conductivity  and  electrode  temperature.  This  value  is  compa¬ 
rable  to  a  similarly  defined  value  based  on  the  Hall  current  density  in 
the  core  of  a  large-scale  channel. 


Figure  3.3  Electrical  conf iguration  used  for  the  main  experiments. 


Laser  Doppler  anemometry  system 

Plasma  velocities  are  measured  using  laser  Doppler  anemometry.  The 
LDA  system,  shown  schematically  in  Fig.  3.4,  uses  a  dual-beam,  single¬ 
color,  backscatter  configuration.  Velocity  measurements  are  made  at 
numerous  locations  over  the  cross-plane  defined  by  the  port  location. 
Translation  of  the  measurement  volume  in  the  z-direction  is  achieved  by 
translating  the  focusing  lens  along  the  optical  axis,  and  in  the  y- 
direction  by  rotating  about  its  horizontal  axis  a  large  mirror  which 
directs  the  converging  beams  toward  the  small  mirror  next  to  the  optical 
port . 

Conductivity  measurements 

The  effect  of  secondary  flow  on  the  plasma  conductivity  profile  is 
investigated  by  driving  a  Faraday  discharge  across  a  set  of  electrodes 
whose  external  circuitry  is  separate  from  that  used  to  drive  the  Hall 
current.  This  set  consists  of  three  electrode  pairs  located  just 
upstream  of  the  optical  port.  These  electrodes  are  used  only  for  the 
conductivity  measurements,  and  are  open-circuited  at  all  other  times. 
The  outer  two  pairs  serve  as  guard  eLectrodes,  so  that  the  Faraday 
current  density  Jy  flowing  across  the  gap  of  the  middle  pair  can 
be  estimated  more  accurately.  Near-electrode  voltage  drops  in  both  the 
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Figure  3.4  Schematic  of  laser  Doppler  anemometry  optical  train. 


x-  and  y-directions  are  obtained  from  a  grid  of  iridium  voltage  pins 
inserted  through  a  channel  sidewall.  Voltage  distributions  are  obtained 
over  a  range  of  values  for  Jy,  and  under  various  applied  conditions, 
including : 


1.  applied  Faraday  current  only; 

2.  both  Faraday  current  and  Hall  current; 

3.  Faraday  current  with  magnetic  field;  and 


4. 


Faraday  current,  Hall  current  and  magnetic  field. 


For  cases  (3)  and  (4)  measurements  are  made  with  the  magnetic  field 
pointing  in  both  the  positive-  and  negative-z  directions.  The  open- 
circuit  voltage  with  magnetic  field  is  measured  as  well.  Note  that  only 
case  (4)  would  be  expected  to  drive  secondary  flow. 

The  voltage  information  is  analyzed  using  an  approximate  model  based 
on  the  expression,  obtainable  from  the  generalized  Ohm's  law. 


E  -  uB  +  BE 


(l  +  b2)j 


(3.4) 


where  p  is  the  resistivity,  and  Ex  the  electric  field  components  and 
u  the  axial  velocity.  The  overbar  indicates  an  average  over  a  region  of 
the  cross-plane  which  encompasses  the  width  W  in  the  z-directlon  and 
ranges  in  the  y-direction  from  the  electrode  surface  to  the  height  of 
the  nearest  voltage  pin,  a  distance  of  6.2  mm,  or  H/8,  H  being  the 
channel  height. 


Temperature  measurements 

All  electrodes  are  fitted  with  thermocouples  located  close  to  their 
surfaces.  Of  particular  interest  are  the  temperatures  of  the  electrode 
pair  just  upstream  of  the  velocity  measurement  plane,  for  cases  in  which 
these  electrodes  are  not  passing  current.  A  number  of  measurements  are 
recorded  for  runs  in  which  a  case  without  secondary  flow  (i.e.,  without 
the  simultaneous  presence  of  magnetic  field  and  Hall  current)  is 
immediately  followed  or  preceded  by  a  case  with  secondary  flow  (i.e., 
both  were  present).  The  change  in  surface  temperature  between  these  two 
cases  can  then  be  attributed  to  altered  heat  transfer  as  a  result  of 
secondary  flow.  As  with  the  other  measurements,  data  is  collected  for 
both  the  positive-  and  negative-z  orientations  ot  the  magnetic 
induction. 


Experimental  results 

Transverse  velocity  measurements 

Non-MHD  case 

So  as  to  have  a  base  case  with  which  to  compare  the  measurements 
with  MHD  effects  present,  measurements  of  the  y-directed  velocity  have 
been  made  for  the  case  of  zero  magnetic  induction  and  zero  current.  The 
results  are  shown  in  Fig.  3.5.  The  peak  measured  transverse  velocity  is 
0.05  U^.  This  is  significantly  higher  than  reported  measurements  for 
secondary  flows  driven  purely  by  the  interaction  of  turbulence  stresses 
with  the  corners  of  a  rectangular  duct.  Several  factors  may  be  involved 
in  this  discrepancy,  including  combustion  non-uniformities  and  other 
non-idealities  related  to  the  combustor,  density  gradients,  and  the 
heterogeneous  nature  of  the  channel  walls  (the  top  and  bottom  walls  had 
alternating  sections  of  brick  and  of  electrodes,  whereas  the  sidewalls 
had  only  brick). 

MHD  cases 


Measurements  of  the  y-directed  velocity  have  been  made  for  a  range 
of  conditions,  and  with  the  magnetic  field  pointing  in  both  the  posi¬ 
tive-  and  negative-z  directions.  Reversing  the  magnet  polarity  reverses 
the  sign  of  the  Lorentz  force,  and  so  should  be  approximately  equivalent 
to  rotating  the  channel  by  180°  about  its  axis. 

The  measurements  for  the  case  of  B  pointing  in  the  positive-z  direc¬ 
tion  are  shown  in  Fig.  3.6.  In  this  figure  and  in  the  figures  following 
the  secondary  flow  appears  stronger  near  the  left  sidewall  than  near  the 
right;  we  believe  that  this  is  an  artifact  caused  by  a  consistently 
erroneous  measurement  of  the  wall  location.  Shifting  all  the  measure¬ 
ment  locations  to  the  left  by  an  amount  A(z/W)  =  0.035  would  produce  a 
symmetric  result,  which  is  physically  what  we  would  expect.  This 
systematic  error  is  additional  to  the  uncertainty  in  spatial  location. 


In  any  case,  the  results  suggest  the  presence  of  two  large,  counter¬ 
rotating  secondary  flow  cells,  with  the  flow  in  the  core  directed  down¬ 
ward  in  agreement  with  the  first-order  argument  illustrated  in  Fig.  3.1. 
The  peak  measured  velocities,  found  in  the  sidewall  region,  measured 
0.12  Ub;  as  the  magnitude  of  the  transverse  velocity  is  increasing 
steeply  toward  the  sidewall,  it  is  reasonable  to  assume  that  signifi¬ 
cantly  higher  velocities  would  be  found  by  measurements  closer  to  the 
walls . 

The  vorticity  is  concentrated  in  the  top  half  of  the  channel,  with 
the  vortex  cells  closing  at  about  y/H  =  0.35.  An  explanation  for  this 
behavior  is  suggested  by  a  three-dimensional  numerical  simulation  per¬ 
formed  by  Maxwell,  Early  and  Demetriades  [3.12]  for  conditions  which 
were  similar  to  those  of  the  present  experiments.  In  these  calculations 
the  secondary  flow  field  at  the  upstream  end  of  the  interaction  region 
has  a  symmetric  structure  similar  to  that  in  Fig.  3.1.  As  the  flow 
progresses,  hot  plasma  from  the  core  has  time  to  be  swept  sideways  and 
then  upwards  as  it  follows  a  secondary  flow  path.  This  has  the  effect  of 
increasing  3Jx/3z  in  the  upper  sidewall  regions,  which  in  turn 
intensifies  the  vorticity  in  these  regions.  The  coupling  between  3Jx/3z 
and  axial  vorticity  causes  the  secondary  flow  cells  to  migrate  toward 
the  top  walls  in  the  downstream  half  of  the  channel.  In  this  respect  the 
simulation  and  the  measurements  are  in  good  agreement.  The  quantitative 
agreement  is  also  reasonably  good,  although  precise  comparisons  would 
require  a  more  detailed  specification  than  is  available  of  the  sidewall 
temperature  distributions  and  of  the  upstream  turbulence  history.  Both 
of  these  factors  were  shown  in  the  simulation  to  have  a  substantial 
influence  on  the  secondary  flow  field:  the  sidewall  temperature  because 
it  affects  SJ^/Sz,  the  turbulence  intensity  because  of  its  cross-plane 
smoothing  effect. 

The  measurements  in  the  region  at  the  bottom  of  the  cross-plane 
indicate  a  flow  which  is  everywhere  directed  away  from  the  bottom  wall. 
A  possible  explanation  for  these  measurements  is  discussed  in  connection 
with  the  temperature  results. 


In  Fig.  3.7  the  results  for  the  case  with  the  magnet  polarity 
reversed  are  shown.  As  anticipated  the  velocity  field  is  upside-down 
compared  to  the  results  in  Fig.  3.6.  A  few  measurements  have  been  made 
in  this  case  closer  to  the  sidewalls,  resulting  in  a  peak  transverse 
velocity  measurement  of  0.15  U^* 

Effects  of  Secondary  Flow 

The  measured  secondary  flow  field  is  certainly  strong  enough  to  have 
an  effect  on  several  aspects  of  the  plasma  momentum  and  thermal 
behavior.  This  section  reports  the  results  of  the  measurements  of  mean 
axial  velocity,  turbulence  intensity,  electrical  conductivity  and  elec¬ 
trode  surface  temperature. 
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Figure  3.7  Measurements  of  the  mean  y-directed  velocity  with 

the  magnetic  field  point  in  the  negative  z-direction 


Mean  axial  velocity 

Measurements  of  the  mean  axial  velocity  with  the  magnetic  field 
pointing  in  the  positive-z  direction  are  shown  in  Fig.  3.3;  these 
results  correspond  to  the  secondary  flow  field  shown  in  Fig.  3.6.  The 
distortion  of  the  velocity  profile  is  dramatic:  the  peak  velocity  is 

forced  downward  from  the  center  of  the  cross-plane  to  a  position  of 
about  y/H  =  0.3.  The  isovels  near  the  top  wall  follow  the  vortex 

pattern,  as  slow-moving  fluid  is  swept  upward  along  the  sidewall  and 
then  pushed  downward  along  the  z-center line .  The  axial  velocity  profiles 
along  the  z-centerline  for  three  conditions  of  the  magnetic  field  are 
shown  in  Fig.  3.9.  As  expected  the  profile  is  symmetric  when  Lorentz 
forces  were  absent;  with  Lorentz  forces  present  the  profile  became 
strongly  skewed  according  to  the  direction  of  the  force. 

These  results  are  similar  to  the  measurements  reported  by  McClaine, 
Swallom  and  Kessler  [3.14],  In  their  experiments,  performed  in  the 
larger  Avco  Mk  VI  channel,  the  conditions  included  a  magnetic  induction 
of  4  T,  an  interaction  length  of  2.5  m,  and  peak  velocities  of  about 
1000  m/s.  Unfortunately  the  electrical  loading  which  is  used  during 
these  particular  measurements  is  not  reported;  their  tests  encompassed  a 
wide  range  of  loading  conditions  in  both  Faraday  and  diagonal  configura¬ 
tions  . 

Turbulence  intensity 

Previous  measurements  of  turbulence  intensity  in  the  Stanford  M-2 
channel  have  found  relatively  high  values  of  the  center-point  axial 
turbulence  intensity  u^/U^,  where  uc  is  the  turbulence  velocity  (the 
standard  deviation  of  the  velocity  probability  distribution  function) 
and  uc  is  the  center-point  mean  velocity.  For  example,  Reis  and  Kruger 
[3.16]  measured  values  ranging  from  6.6%  to  9.b %  for  a  variety  of  non- 
MHD  flow  cases;  in  the  present  experiments  several  measurements 
of  u  /u  yielded  values  ranging  from  5.8%  to  6.3%.  For  comparison, 
Melllng  and  Whitelaw  [3.15]  in  their  water  flow  experiment  measured 


values  of  4-5%  ac  an  axial  location  where  the  boundary  layers  had 
merged.  It  has  been  conjectured  that  the  relatively  high  axial  turbu¬ 
lence  in  the  Stanford  channel  is  caused  by  combustion  nonuni forrai t ies . 
The  measured  value  in  the  non-MHD  case  of  the  y-directed  core  turbu¬ 
lence  v,  again  normalized  by  uc,  is  3.75%;  this  is  virtually  identical 
to  the  value  reported  by  Melling  and  Whitelaw  [3.15]. 

Measurements  of  the  y-directed  turbulence  intensity  in  the  non-MHD 
case  are  shown  in  Fig.  3.10,  in  which  the  measured  values  of  v  have  been 
normalized  by  U^.  As  expected  the  distribution  is  symmetric  and  |v|  is 
greater  near  walls  parallel  to  y  than  near  walls  normal  to  y. 

The  measurements  made  in  the  MHD  case,  with  the  magnetic  field 
pointing  in  the  positive-z  direction,  are  shown  in  Fig.  3.11.  The  region 
of  minimum  turbulence  (whose  magnitude  is  somewhat  less  than  in  the  non- 
MHD  case  because  of  turbulence  damping  by  the  magnetic  field;  cf.Reis 
and  Kruger  [3.16])  is  shifted  downward  by  Lorentz  forces. 

The  measurements  for  u/U^,  shown  in  Fig.  3.12,  exhibited  the  same 
effect;  these  same  measurements  are  shown  in  Fig.  3.13  normalized 
instead  by  the  local  mean  velocity  u.  In  these  measurements  there 
appears  to  be  a  region  beneath  the  center  of  the  top  wall  where  the 
turbulence  is  greater  than  its  value  closer  to  the  wall;  this  might  be 
explainable  in  terms  of  secondary  flow  convection. 


Electrical  conductivity 


The  results  of  the  conductivity  measurements  for  the  MHD  cases  are 
summarized  in  Fig.  3.14,  in  which  the  near-electrode  resistivity  defined 
by  equation  (b)  has  been  normalized  by  the  value  of  the  core  conduc¬ 
tivity  measured  without  magnetic  field.  The  terms  "anode  region"  and 
"cathode  region"  refer  respectively  to  the  regions  adjacent  to  the  top 
and  bottom  electrodes  which  were  used  for  the  conductivity  measurements. 
The  abscissa  in  this  figure  represents  the  value  of  the  core  current 
density  associated  with  the  Faraday  discharge;  for  all  the  data  shown 
the  Hall  circuit  is  active  as  well,  corresponding  to  case  (4)  described 
earlier . 
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Figure  3.14  Summary  of  tesistivity  results  for  cases  with  complete 
electrical  configuration  and  magnetic  field. 

The  decrease  in  resistivity  as  Jy  increased  is  an  expected  conse¬ 
quence  of  Joule  heating.  Otherwise  these  results  indicate  that  secondary 
flow  caused  the  resistivity  distribution  to  become  significantly  asym¬ 
metric.  The  boundary  layer  receiving  cold  sidewall  fluid  is  much  more 
resistive  than  the  boundary  layer  toward  which  hot  core  fluid  is 
directed;  the  average  value  of  the  ratio  of  the  higher  resistivity  to 
the  value  at  the  opposite  wall  is  1.8.  The  importance  of  these  results 
is  that  near-electrode  voltage  drops,  which  are  a  performance-limiting 
effect  in  combustion  MHD  devices,  can  be  significantly  altered  by  sec¬ 
ondary  flow. 

Electrode  Surface  Temperatures 

Thirteen  separate  measurements  have  been  made  of  the  effect  of 
secondary  flow  on  electrode  surface  temperatures,  as  described  earlier. 
The  results  are  shown  in  Fig  3.15.  Each  point  in  this  figure  represents 
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Figure  3.15  Measured  electrode  surface  temperatures  with 

Hall  current  and  magnetic  field  minus  reference 
measurements  without  current. 


the  temperature  with  secondary  flow  (with  simultaneous  Hall  current  and 
magnetic  field)  minus  the  temperature  without  secondary  flow  (without 
current);  the  latter  temperature  is  typically  about  1000  K.  A  separate 
analysis  of  data  not  shown  in  the  figure  indicates  that  the  presence  of 
magnetic  field  alone  makes  no  difference  in  the  temperature,  while  the 
presence  of  Hall  current  alone  causes  a  warming  of  both  the  top  and 
bottom  electrodes  by  an  average  of  about  10  K,  because  of  Joule  heating. 

The  effect  of  secondary  flow  in  the  positive-B  case  is  seen  to  be  a 
substantial  heating  of  the  bottom  electrode  and  a  smaller  but  still 
significant  cooling  of  the  top  electrode;  reversing  the  magnet  polarity 
reversed  the  results.  The  electrode  toward  which  hot  core  fluid  is 

pushed  is  heated  on  the  average  by  85  K;  the  opposite  electrode  is 

cooled  by  an  average  of  32  K.  Using  the  fact  that  the  electrodes  are 

cooled  by  water  at  about  300  K,  we  calculate  that  secondary  flow  evi¬ 
dently  causes  a  12%  increase  in  the  heat  flux  from  the  plasma  to  the 
heated  electrode  and  a  5%  reduction  in  the  heat  flux  to  the  cooled 

electrode. 
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Aside  from  the  obvious  conclusion  that  secondary  flow  has  Che 
expected  qualitative  asymmetric  effect  on  heat  transfer  to  the  top  and 
bottom  walls,  it  is  interesting  to  consider  the  surprising  result  that 
the  effect  at  the  heated  electrode  is  more  than  twice  as  great  as  the 
effect  at  the  cooled  electrode.  This  result  is  surprising,  first, 
because  the  secondary  flow  cells  appear  to  be  concentrated  near  the 
cooled  electrode,  and  second,  because  the  resistivity  measurements  show 
the  opposite  trend,  with  a  stronger  effect  being  observed  in  the  region 
adjacent  to  the  cooled  electrode. 

A  possible  explanation  is  as  follows.  As  already  noted,  the  second¬ 
ary  flow  cells  not  only  sweep  cold  sidewall  fluid  into  the  top  wall 
boundary  layer,  but  also  eventually  convect  hot  core  fluid  around  a  loop 
into  the  same  boundary  layer  in  the  sidewall  region.  As  the  stainless 
steel  electrodes  are  thermally  quite  conductive,  the  net  result  is  that 
the  cooling  effect  is  mitigated;  no  such  effect  occurs  for  the  boundary 
layer  of  the  bottom  wall,  toward  which  core  fluid  has  been  pushed  but 
which  has  less  vortical  mixing.  Distortion  of  the  turbulence  distribu¬ 
tion  could  also  be  a  factor  in  the  asymmetrical  heat  transfer.  Near¬ 
electrode  voltage  drops,  on  the  other  hand,  are  dependent  on  complex 
physical  processes  which  are  strongly  but  non-linearly  temperature- 
dependent.  Since  the  calculated  resistivities  have  been  averaged  over  an 
area  described  earlier,  it  is  not  remarkable  that  the  net  quantitative 
effects  for  the  resistivity  distribution  could  be  quite  different  than 
for  heat  transfer,  although  the  qualitative  trends  are  the  same. 

The  temperature  results  suggest  a  possible  explanation  for  the 
secondary  flow  measurements  in  the  region  near  the  bottom  wall  (in  the 
positive-B  case).  The  substantial  warming  of  the  bottom  electrode 
suggests  that  the  mass  density  near  the  bottom  wall  may  have  been 
decreasing  in  the  axial  direction,  causing  a  net  reduction  in  the  axial 
mass  flux  which  would  have  to  be  balanced  by  a  flow  away  from  the  wall. 
Although  this  behavior  did  not  appear  in  the  simulation  of  Maxwell, 
Early  and  Deraetriades  [3.12],  their  calculations  assumed  that  the  wall 
temperature  is  constant. 


Conclusions 


In  conclusion,  che  past  year  has  been  one  of  significant  achievement 
for  our  program  in  MHD-induced  secondary  flow.  The  results  constitute 
the  first  direct,  quantitative  measurements  of  secondary  flow  in  a 
combustion  MHD  channel.  Our  experiments  are  showing  that  strong 
secondary  flows  can  be  driven  by  Lorentz  forces  which  are  present  in 
moderate-to  high-interaction  MHD  devices,  and  that  these  secondary  flows 
can  have  a  profound  impact  on  the  performance  of  these  devices. 
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